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Abstract 

The scaling properties of self-avoiding polymerized 2-dimensional membranes are 
studied via renormalization group methods based on a multilocal operator product 
expansion. The renormalization group functions are calculated to second order. 
This yields the scaling exponent v to order . Our extrapolations for v agree with 
the Gaussian variational estimate for large space dimension d and are close to the 
Flory estimate for d = 3. The interplay between self- avoidance and rigidity at small 
d is briefly discussed. 
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1 Introduction 



The statistical properties of polymerized flexible membranes are an interesting subject [1], 
which is still only partially understood. These objects, also called tethered membranes, 
are two-dimensional generalizations of polymers (one-dimensional flexible chains). It is 
expected that such two-dimensional membranes exhibit a larger variety of behavior than 
polymers, when the temperature and the elastic properties of the membranes are varied. 
One reason is the following: simple dimensional analysis shows that for two-dimensional 
films, the bending rigidity modulus has the dimension of a pure energy and is there- 
fore marginally relevant (in the sense of the renormalization group) [2]. Moreover, two- 
dimensional membranes may have internal shear elasticity and this separates two very 
different classes of flexible membranes: fluid membranes with zero shear modulus and 
crystalline (or tethered) membranes with non-zero shear modulus [3]. 

In this paper, we consider tethered membranes. As long as one takes into account 
only the local elastic forces (i.e. bending rigidity, compression and shear moduli), numer- 
ical simulations and analytical calculations point towards a consistent and relatively well 
understood picture [4, 5, 42, 6, 7, 8]. For high rigidity or equivalently low temperature, 
the membrane is in a flat phase with an average orientation and a classical "fractal" di- 
mension dp = 2. The roughness properties of this flat phase are nevertheless non-trivial, 
due to the non-linear coupling between undulation modes and phonons. For low rigidity 
or equivalently high temperature, the membrane is in a crumpled phase, without any 
global orientation. In this phase the statistics of the surface is Gaussian (at large length 
scales) and its fractal dimension is inflnite {dp = oo). The two phases are separated by a 
crumpling transition, which occurs for a flnite bending rigidity (or equivalently a non-zero 
temperature). Numerical simulations and some of the analytical calculations indicate that 
this transition is continuous and characterized by non-trivial critical exponents. 

The above results do not take into account steric interactions, that is local self- 
avoidance and thus concern "phantom surfaces" . For membranes {D = 2) these interac- 
tions must be relevant for the crumpled phase (at large length scales), whatever the bulk 
dimension d of space is, since the fractal dimension of phantom surfaces is then infinite. 
They are expected to be strongly relevant for physical membranes in three-dimensional 
space {d = 3). This is in contrast with polymers, where steric intercations are relevant 
only for d < A and lead to a swollen phase in dimensions 1 < d < 4 with 1 < dp < 2. 

The study of the effect of self-avoidance for tethered membranes is in fact much more 
difficult than for polymers. Most of the studies rely on numerical simulations, which find 
evidence for a flat phase in 3-dimensions [9, 19]. Some experiments have been performed 
using thin sheets of graphite oxide, but the results are contradictory. In [24] but not in 
[25] a swollen crumpled phase is found. 

In order to study theoretically self-avoiding membranes, one can use standard approx- 
imations, such as the Flory approach or variational methods. A more systematic renor- 
malization group approach has been initiated in [27, 28]. This approach is inspired from 
the direct renormalization method used for polymers [29]. It is a perturbative method 
and it has been used to calculate the fractal exponent u* = 2/dp at flrst order in an e- 
expansion [27, 28, 30]. An important problem has been to check the internal consistency 
of this renormalization method and to extend it properly to all orders in perturbation 
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theory. This has been completed by B. Duplantier, E. Guitter and one of the authors 
in [31] (a detailed proof might be released soon [32]), who have shown that the model of 
self-avoiding membranes of [27, 28], although corresponding to a non- /oca/ field theory, is 
renormalizable in perturbation theory. This result establishes the validity of the 0{e) cal- 
culations and gives a systematic formalism to extend these calculations to higher orders. 

In this paper we present for the first time in full details the renormalization group 
calculation at second order for the model of self-avoiding tethered membranes and discuss 
the results obtained by this approach. A short presentation of the main results has already 
appeared in [33]. 

This study is interesting for several reasons: 

(1) It provides an explicit realization of the renormalization formalism of [31] . In particular 
these calculations illustrate nicely how subdominant divergences are organized. 

(2) The calculation is a non-trivial task. As usual in renormalization theory, there is a 
long way from existence and convergence theorems to actual calculations, but additional 
difficulties arc present here. In particular, beyond first order, the amplitudes (which are 
called manifold integrals and generalize in a non-trivial way Feynman integrals in the 
Schwinger representation to non-integer dimensions < D < 2) cannot be calculated 
analytically and in fact describe distributions. One must rely on numerical estimates for 
these integrals. 

(3) It leads to the first estimates at order 0{e'^) for the fractal exponent u* = 2/ dp. 
This allows to check the validity of the 0{e) estimates, as well as the consistency of the 
extrapolation methods used to extract these estimates. Indeed, direct calculations for 
two-dimensional membranes are impossible, but the model is extended to "membranes" 
with internal dimension < D < 2, thus interpolating between polymers {D = 1) and 
membranes {D = 2). Perturbative calculations lead to an £-expansion, where the e- 
parameter is given by 

£ = £{D,d) = 2D-d^^^—^. (1.1) 

D is the internal dimension of the membrane and d the dimension of bulk space (in 
which the membrane fiuctuates). Using the e-expansion one may start a-priori from any 
point (Do,c?o) such that ^(Doi'^o) = to extrapolate to (for instance) the physical point 
(Z) = 2, d = 3). In [30] an extrapolation scheme was used for 0{e) calculations. We shall 
need and will develop more systematic extrapolation schemes for 0{e^) calculations and 
we shall discuss their respective advantages. 

(4) Finally, the calculations can be performed for membranes {D = 2) in a space with 
arbitrary dimension d and we can compare explicitly our 0{e^) results for v* with other 
predictions. It turns out that our estimates are quite reliable for large d and in remarkable 
agreement with the result of a Gaussian variational estimate. We shall explain this fact 
and argue that this feature persists at higher order in perturbation theory. For smaller d 
the results are less stable, but still good and in reasonable agreement with Flory estimates 
(we have no good explanation for this fact). 

The paper is organized as follows: 
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In section 2 the continuous model of self-avoiding tethered membranes of [27, 28] is 
introduced. It is a generahzation of the Edwards-model for polymers. We recall the basic 
results of [31, 32] concerning the structure of the perturbative expansion for the model, 
the nature of the short distance, or ultra-violet (UV), divergences and their relation with 
the so-called multilocal operator product expansion (MOPE). 

In section 3 we recall the renormalization of the model at first order in perturbation 
theory (one loop), give the explicit expressions for the counterterms in the "minimal 
subtraction scheme" used in this paper and derive the renormalization group functions 
(3{b) and z/(6) at one loop. The results are of course not new, but this fixes the method 
and the notations used throughout the paper. 

In section 4 we analyze the UV-divergences at second order (two loops), and obtain 
the expressions for the two-loop counterterms, in terms of singular parts of integrals of 
MOPE coefficients. First the UV-singular configurations (divergent and subdivergent dia- 
grams) are identified (subsection 4.1), then the leading UV-divergences (poles in 1/e^) are 
obtained and their exponentiation (predicted by renormalization group) checked (subsec- 
tion 4.2). This allows to obtain the subleading divergences (poles in l/e), whose residues 
give the renormalization group functions at two loops (subsection 4.3). We then express 
these residues as combinations of convergent integrals involving MOPE coefficients (sub- 
section 4.4). 

The next seven sections are devoted to the explicit calculation of these integrals. 
In section 5 we briefly introduce some basic analytical methods used in the calculations, 
which have been developed in [34]. These are: (1) the definition of the "distance measure" 
which allows to define properly the integration over non-integer dimensional space IR^ 
{0 < D < 2); (2) the expression of the residue of the UV-poles in {l/s) as a boundary 
term in the distance integrals; (3) the "conformal mapping" technique, which allows to 
map different domains of integration and leads to crucial simplifications. 

Sections 6, 7 and 8 are devoted to the numerical calculation of the three diagrams 
which contribute to the coupling constant renormalization. In all cases, we start from 
the corresponding MOPE cocfiicicnts and determine explicitly the integrals which have 
to be computed (this is in general not straightforward). Then we evaluate numerically 
these integrals for values of 1 < D < 2 (the internal dimension of the membrane). We 
have in general to decompose the domain of distance integration into several pieces, called 
sectors, and to find ad hoc changes of variables in each sector. In addition, this requires 
an adaptive Monte Carlo integration routine, first developed in [34], in order to master 
the rapid variations of the integrand. For the diagram of section 7 there are additional 
subtleties, arizing from the fact that the measure of integration is then a distribution, with 
non-integrable singularities on some boundary of the integration domain, which have to 
be treated by a finite part prescription. For each diagram we compute analytically its 
D — > 1 limit. This provides a check of the numerics. Secondly for D = 1 our model 
reduces to the Edwards-model for polymers, for which 2-loop calculations have already 
been performed by several authors and which give an additional check of our calculations. 

Sections 9, 10 and 11 are devoted to the numerical calculation of the three diagrams 
which contribute to the field renormalization and are organized in a similar way. 

In section 12 we use these two loop results to calculate critical exponents for self- 
avoiding membranes. First we recall how the renormalization group functions at two 
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loops are related to the counterterms that we have calculated (subsection 12.1) and thus 
obtain the £^ term for the fractal exponent v*. Then as explained above, we have to set up 
extrapolation methods to extrapolate from the e-expansion to the physical case D = 2. We 
generalize and systematize the extrapolation method proposed by Hwa [30], and use our 
new schemes to evaluate i^* for membrane for various bulk space dimensions 2 < d < oo 
(subsection 12.2). In subsection 12.3 our results are compared with that of the Gaussian 
variational method. We argue that u*, as obtained from the (properly resummed) e- 
expansion, must coincide with the variational estimate z/yar for large d and we propose a 
new e-expansion for u* which coincides with z/var at order and which is shown to give 
very good results for large d. In subsection 12.4 we compare our results for u* with that 
of the Flory method and consider a similar new £-expansion around z^Fiory The results 
of all the two loop extrapolations for i/* are summarized in subsection 12.5. Finally we 
briefly discuss the case of other scaling exponents for self-avoiding membranes, namely the 
correction to scaling exponent u and the contact exponent 62 (subsection 12.6), as well as 
the fractal exponent i^^ for membranes at the tri-critical ^-point, which has been already 
calculated at order e. In subsection 12.8 we briefly summarize the main results from 
numerical simulations as well as experimental data. We also present a heuristic argument 
which explains why both in numerical simulations and experiments self-avoiding tethered 
membranes are found in a fiat phase (subsection 12.9). 

Conclusions and future prospects are given in section 13. Several technical points or 
examples are gathered in the appendices. 

2 Definition of the Model 

We start from the continuous model for a D-dimensional flexible polymerized membrane 
introduced in [27, 28]. This model is a simple extension of the well known Edwards' model 
for continuous chains. The membrane fluctuates in rf-dimensional space. Points in the 
membrane are labeled by coordinates x G IR^ and the conflguration of the membrane in 
physical space is described by the field f : x e — > f{x) e IR'^. The free energy for a 
configuration is given by the bare Hamiltonian 

n[r\ = ^ ^(Vr»)^ + bJJ^ ~S'{f{x) - f{y)) . (2.1) 

The integral runs over D-dimensional space and V is the usual gradient operator. The 
normalizations (hidden in -^^d-i Jx 5^{r — f) ) are chosen in order to simphfy the 
calculations, but are unimportant for the general understanding (see appendix A). The 
first term is a Gaussian elastic energy which is known to describe the free "phantom" 
surface. The interaction term corresponds to a weak repulsive contact interaction (for 
6 > 0). The expectation value of physical observables are obtained by performing the 
average over all field-configurations f{x) with the Bolzmann weight e"^'^. 

Perturbation theory is constructed by performing the series expansion in powers of the 
coupling constant h. This expansion suffers from ultraviolet (UV) divergencies which have 
to be removed by renormalization and which are treated by dimensional regularization, i.e. 
analytical continuation in D and d. A physical UV-cutoff could be introduced instead but 
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would render the calculations more complicated. Long-range infrared (IR) divergencies 
also appear. They can be cured by using a finite membrane, or by studying translationally 
invariant observables, whose perturbative expansion is also IR-finitc in the thermodynamic 
limit (infinite membrane). Examples of such observables are "neutral" products of vertex 
operators 

N N 

O ^ H : g'^'^^^''^) : ^4 = 0. (2.2) 

a=l a=l 

In the following we discuss the renormahzation of the model, i.e. we only deal with the 
UV-divergencies. 

Let us first analyze the theory by power-counting. In internal momentum units, such 
that [x] — —1, the dimension of the field and of the coupling-constant are: 

2 — D 

[r] = -v = — , \^\=e = '2D -vd (2.3) 

In the sense of Wilson the interaction is relevant for £ > 0. Perturbation theory is then 
expected to be UV-finite except for subtractions associated to relevant operators. We 
will come back to this point later. For clarity we shall represent graphically the different 
interaction terms which have to be considered. The local operators are 

1 = . (2.4) 
\{^T{x)f = ^ . (2.5) 

The bi-local operator, the dipole, is 

S'^irix) -r{y)) = • .. (2.6) 

The expectation- value of an observable is: 

fV\r] 0\r] e~'^M 

Perturbatively, all expectation-values are taken with respect to the free theory: 

JV[r]e 2-dJx2\^^) 
A typical term in the expansion of (2.7) is 



where the integral runs over the positions of all dipole-endpoints. The analysis in 
shows that UV-divergencies appear when some dipole-endpoints approach each other. 
The divergencies are analyzed via a multilocal operator product expansion (MOPE). 
The principle is examplified in figure 2.1. One considers n dipoles (here n — and one 
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Figure 2.1: Example of a contraction 



separates the 2n end-points into m separate subsets (here m = 3) dehmited by the dashed 
hnes. The MOPE describes how the product of these n dipoles behaves when the points 
inside each of the m subsets are contracted towards a single point z. The result is a 
sum over multilocal operators . . . , -z^}, depending on the ui points Zi, ■ ■ ■ , z^n, of 

the form 

J2C^'^Ul:Z2,...,Zm} (2.10) 

where the MOPE-coefficients C$ depend only on the relative distances Xi — Xj between the 
dipole end-point positions Xi inside each subset. This expansion is valid as an operator- 
identity, i.e. inserted in any expectation value and in the limit of small distances between 
contracted points. As the Hamiltonian (2.1) does not contain any mass-scale, the MOPE- 
coefficients are homogenous function of the relative positions between the contracted 
points. The degree of homogeneity is given by simple dimensional analysis. In the case 
considered here, where n dipoles are contracted to an operator this degree is simply 
(for the definition of u see (2.3)) 

degree[C^] = -nud + (2.11) 

where [$] is the canonical dimension of the operator $ and d{2 — D)/2 is simply the 
canonical dimension of the dipole. 

In order to evaluate the associated singularity, one has finally to integrate over all 
relative distances inside each subsets. This gives an additional scale factor with degree 
D{n — m). A singular configuration, such as depicted on figure 2.1 , will be UV-divergent 
if its degree of divergence, defined as D{n — m) + deg[C|J^], is negative. It is superficially 
divergent if the degree is zero and convergent otherwise. 

The power counting analysis of [31, 32] shows that at the critical dimension £ = the 
identity operator • is relevant, while the local operator -j- and the bi-local dipole operator 

• • are marginally relevant. Contractions of n dipoles to • give relevant divergences 

(negative powers of the short distance cut-off, or poles for some £ > 0), while contractions 

towards -\- and • • give superficial divergences (logarithms of the cut-off, or poles at 

e = 0). 

The analysis of the renormalization in [31, 32] shows that in order to make perturbation 
theory UV-finite at e = 0, one has to add to the Hamiltonian counterterms proportional 
to all relevant and marginal operators. These counterterms are obtained from the MOPE- 
coefficients, integrated with some IR-cutoff procedure. 
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3 Renormalization at 1-loop order 



Let us continue on the concrete example of the one-loop divergences, where we shall also 
fixe our notations for the MOPE-coefficients. When the end-points {x, y) of a single dipole 
are contracted to a point (taken here to be the center-off- mass z = {x + y)/2), the MOPE 
is 



' y — yx y ' I* 

with the first MOPE-coefficients given explicitly by 



a/3 a/3 + 



(3.1) 



\x-y\ 



-ud 



-^c.p)=-\\x-yr^''^'\x-y)a{x-y)p (3.2) 



and where -j- ap denotes the local tensor operator 



a/3 = ^dafdpr 



The integral over the relative distance x — y for 



(3.3) 



-j- is at £ = logarithmically 



divergent. 

The simplest contraction to a dipole is when two dipoles collapse. The corresponding 
MOPE-coefficient is 



-d/2 



(3.4) 



where x and y are now the relative distances inside the two subsets. Another possibility 
is to consider the contraction \ / . But as 



(3.5) 



X y z ■ • yi ~y'z 

this does not give a pole at £ = 0, but a term proportional to the relevant divergence 
(that we discuss now), times a regular contribution. 

In the next step counterterms have to be introduced in order to subtract these diver- 
gencies. We have to distinguish between the counterterm for the relevant operator and 
those for marginal operators. The first one can be defined by analytic continuation, the 
latter require a subtraction scale. Indeed, the divergence for • is given by the integral 



K-^<\x-y\<L 



^ dx 



-X 



D-i/d 



A-i X 



D 



(3.6) 



where A is a high-momentum UV-regulator and L a large distance regulator. For e ~ 
this is UV-divergent but IR-convergent. The simplest way to subtract this divergence is 
therefore to replace the dipole operator by 



(3.7) 
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where • • = \x 



y\ This amounts to add to the bare Hamiltonian (2.1) the 



UV-divergent counterterm 



y 



—vd 



(3.8) 



X Jy 



which is a pure number and thus does not change the expectation-value of any physical 
observable. This prescription is sufficient to subtract all relevant UV-divergences in the 
calculation of the renormalization-group functions at 2-loop order, that we present in the 

next sections. 

We now treat marginal operators. Let us come back to the MOPE (3.1). The integral 



over the relative distance of 



/ {A 

Jx—y \ < 



-j- a^s ) a^a is logarithmically divergent at e 



0. In order to find the appropriate counterterm, we use dimensional regularization, i.e. set 
£ > 0. An IR-cutoff L, or equivalently a subtraction momentum scale ^ = L~^, has to be 
introduced in order to define the subtraction operation. As a general rule, let us integrate 
over all distances appearing in the MOPE-coefficient, bounded by the subtraction scale 
L — This projects the tensor operator -j-a/j onto the scalar times the integral 



x—y\ <L 



H^)^ = L^f(s,D) 



(3.9) 



Following [31, 32, 34] we use a minimal subtraction scheme (MS). The internal dimension 
of the membrane D is kept fixed and (3.9) is expanded as a Laurent series in s, which 
here starts at e~^. Denoting by ( | the term of order of the Laurent expansion of 
( I )^ for L = 1, the pure pole part of (3.9) is found to be 



1 1 
'2De 



(3.10) 



It is this pole-term that we are going to subtract in the MS-scheme. This is done by 
adding to the Hamiltonian a counterterm 



An[r] 



(3.11) 



Similarly, the divergence arising from the contraction of two dipoles into a single dipole 
is subtracted by a counterterm proportional to the single pole of 



x\<L J\y\<L 




>] = [ I (kr+bp^)" 

) J\x\<LJ\y\<L^ ' 



-d/2 



(3.12) 



As a result the model is UV-finite if we use the renormalized Hamiltonian Tin 

-HniA = ^ I ^(Vr»)^ + hZ,^^ JJ/in^) ' %)) , (3-13) 

instead of the bare Hamiltonian Ti.[r\. Now r and b are the renormalized field and coupling 
constant, and // = is the renormalization momentum scale. The renormalization 
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factors are at one loop 



Z^l-{2-D) 



b + 0{b'') 
b + 0{b^) 



(3.14) 
(3.15) 



The renormalized theory can be reexpressed in terms of the bare (unrenormahzed) 
theory through 

fo(x) = Z^/^r{x), bo = bZbZ'^/^ii' (3.16) 

Following the analysis of [31], the renormalization group /5-function and the anomalous 
scaling dimension v oif are obtained from the variation of the coupling constant and the 
field with respect to the renormalization scale /i, keeping the bare couplings fixed. They 
are written in terms of Z and Z^ as 



d_ 



-eb 



v{b) 



bo 

2-D 



l + 6|lnZ, + f6|lnZ 



_1 d_ 



bo 



(3.17) 
(3.18) 



4 Derivation of the Counterterms at 2-Loop Order 

4.1 The 2-loop counterterms in the MS scheme 

In this section we apply the formalism explained above to determine the counterterms 
which renormalize the theory at second order. If we consider only the bare theory, given 
by the bare Hamiltonian (2.1), power counting gives the three UV-divergent diagrams 
(together with their symmetry factors) 




which give short distance singularities when the points inside the subsets are contracted 
to a single point. These singularities give double and single poles at £ = 0. There are two 
other potentially dangerous diagrams 




These diagrams do not give poles at £ = for reasons similar to what happens with (3.5). 
Now one has to remember that the model is already renormalized at 1-loop, i.e. that we 
use the renormalized Hamiltonian (3.13), with the counterterms (3.14) and (3.15). As a 
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consequence there are five additional divergent diagrams, which come from the insertion 
of the 1-loop counterterms 



1^ 




(e) 




if), -2//- 




(9) 



2^ 



-2£ 




{h) 



(4.3) 



There are other potentially divergent diagrams, analogous to those depicted in (4.2), 
which factorize into convergent diagrams. 

The first four terms in (4.3) arc a combination of a diagram divergent at 1-loop order 
(giving a single pole) times a divergent 1-loop counterterm (which gives another single 
pole). The fifth term is more peculiar: it is the combination of a convergent diagram 
(which corresponds to a contact term) times two 1-loop counterterms (thus giving also a 
double pole). 

Owing to the MOPE, diagrams (a), (e), (/) and {h) give a divergence proportional 
to the insertion of the local operator . With the notations introduced in the previous 
section, they can be subtracted by adding a counterterm proportional to the divergent 
part of the integral of the corresponding MOPE coefficients 






1 



■2e 




(4.4) 



Since we use the minimal subtraction scheme, we want to subtract only the double and 
single poles in e at e = 0. To isolate these poles, we have to perform a Laurent expansion 
of the various terms in (4.4) and to keep the terms of order and but to drop the 
analytic part. Setting the renormalization momentum scale // = L~^, we obtain the final 
expression for the renormalization factor Z at 2-loop order 



Z = 1 - 6(2- D) 
+ 6' (2-^) 



2(0^^0 
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1 

-. + 2 



(4.5) 



Here ( | )^„j denotes the sum of the terms of order s"^, • • • , e'^" in the Laurent expan- 
sion of ( I taken at L = 1. 

Similarly, the diagrams (6), (c), (e) and (/) give a divergence proportional to the 

bilocal operator • • . An analogous analysis leads to the following expression for the 

coupling-constant renormalization-factor at 2 loops 



1 + b 



+ 2 





+ 0{b^) 



(4.6) 



4.2 The leading divergences (double poles) 

In fact we are interested only in the residues of the single poles, that is into the residues 
Ci and /i in the Laurent expansion of the counterterms 



ei 



7i , /: 



l + -b+\^ + ^\b' + 0{b') 



e 



Z,= l + ^(,+ (^ + ^)6^ + 0(6'). 



(4.7) 



Indeed, the finiteness of the renormalization group functions (3.17) and (3.18) implies 
that the residues of the double poles, ji and C2, can be expressed in terms of the 1-loop 
residues e\ and a\. 

Let us show explicitly how this happens, since this will be useful in order to subtract 
efficiently the double poles in the counterterms. From the analysis at 1-loop order, we 
know that 



7 = -(2-i^; 

£ 



Oi 

e 



and we obtain the 1-loop RG functions 

+ ai 



/3(6) = -£6 + 62(-^^i 



2 2 



(4.8) 



(4.9) 
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Inserting these 1-loop RG functions into (3.17) and (3.18), yields differential equations 
for the counterterms Z and Z^, whose solution is explicitly 



1 - 



h ( dci 



d 6]^ +2 



+ ai 



1 - 



h ( dci 



d +2 



+ ai 



(4.10) 



This shows the exponentiation of the leading divergences, i.e. that the residue of the 
leading pole in at order 6" is determined by the residue of the simple pole in at 
order h. Expanding (4.10) to order 6^, we obtain 



1 + ^6 + + 0{h') , 



e e 
with, in terms of diagrams. 



X 



(4.12) 



C2(£) 




— ud 



(4.13) 



Both /2(£) and C2{e) contain not only a constant term but also a term linear in e. The 
latter is not fixed by the RG-functions at 1-loop order. 

The leading poles of Z^^^ (resp. Z^'') must be equal to those of Z (resp. Z},). This 
implies that 

/2 = h{e) + 0{e) ■ C2 = C2(£) + C»(£) (4.14) 

This is compatible with the explicit expressions (4.5) and (4.6) for Z and Z^. Indeed, 
one can directly calculate the double poles in (4.5) and (4.6). The simplest case is the 
diagram (6) in (4.1). A subdivergence occurs when two dipoles are contracted to a single 
dipole. When this contraction is performed first, the MOPE coefficient factorizes as 



(4.15) 




There are three different sub divergences and one finally obtains that the double pole 
associated with this diagram is given by 






(4.16) 



The factor 1/2 comes from the nested integration [35]: the double pole results from the 
integration over a "sector" where the distances inside the subdiagram are smaller than 
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all the other distances. This will become clear in the exphcit calculations of the next 
sections. 

Similarly, let us consider the diagram (c) in (4.1). A subdivergencc occurs when the 
single dipole to the right of the diagram is contracted to a point. The MOPE coefficient 
factorizes as 



(4.17) 




Consequently, the double pole for this diagram is 

_ 1 

-2 ~ 2 





(4.18) 



where the factor 1/2 again comes from the nested integration. 

Finally, let us consider the diagram (a) in (4.1). Four sectors contribute to the double 
pole, which correspond to the subcontractions depicted here 





(4.19) 



Each of the contractions appears with a combinatorial factor two. The double pole for 
this diagram is therefore 




+ 2 X 





(4.20) 



Inserting (4.16), (4.18) and (4.20) in the explicit expressions for the two-loop coun- 
terterms (4.5) and (4.6) we apparently do not obtain for the double poles the results 
(4.12) and (4.13) predicted by the renormalization group. However we can make use of 
the equation of motion to compute the effect of the insertion of the operator -j- in the 
counterterms. In appendix C we show that 

-1 ~ ~Y 



C>(£°) 



(4.21) 




+ 




^+ 0{£') (4.22) 
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Using these identities one recovers (4.12) and (4.13). 



4.3 The subleading divergences (single poles) 



We can now give tiie expressions for the residues of the single poles. For the single pole 
of Z we find 

\ 2 



/i _ (2 - D) 
e 2 



- 2 




^ -2 




One can simplify this expression, since the explicit calculation shows that 



0, 



For the single pole of Z}, the result is 

2 
"3 



Cl 

e 



+ 



-0 



• • 



(4.23) 



(4.24) 



+2 





(4.25) 



4.4 Expressing the residues as convergent integrals 

The next step is to evaluate the residues /i and Ci, that is to write them as convergent 
integrals involving combinations of MOPE coefficients. It is convenient not to compute 
directly /i and Ci, but rather to consider the 0{e^^) parts of the counterterms that are 
not already contained in the second order resummation of the 1-loop divergences, as given 
by (4.12) and (4.13). We thus denote 



C2 _^ Cl 

e 



£^ e 



(4.26) 
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Using d = AD/{2 - D) - 2e/{2 - D) it is easy to extract the 0{e-^) part of (4.12) and 
(4.13). With (4.8) we get 



h = h 



Ci = Ci - 



ei ai 



2(2 -D) ' 2(2 -D) 

Now subtracting (4.12) from (4.5) yields an equivalent exphcit expression for /i 



(4.27) 





v{d + 2) 



(4.28) 



Finally, one has to remember that 



are the terms of order of the Laurent 



series of the integral over distances of the corresponding MOPE coefficient 



L J distances < L 

Using this fact one obtains the following decomposition 

I — U e 

where each term can be written as a convergent integral. These terms are: 



2(G^'0 

_1 / /. 



^\--2 




^2= - 



X 



X 




{2-D){d + 2) 



(4.29) 



(4.30) 



(4.31) 



L 

(4.32) 
(4.33) 
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Even if not written explicitly, we will only calculate the residue of ^i, . . . , ^3 at L = 1. 
Similarly, subtracting (4.13) from (4.6) we obtain for ci 



Cl 

e 




+2 
-2 



that we decompose as 



with 



Co = -2 




2 ^ i2-D) d/, 



) +2 

-2 ^-1 



£-l,£0 3 



Cl 



= Ci+C2+C3 + 0(£°) , 



• + 



(2 - D)d / ; 



(4.34) 
(4.35) 

(4.36) 
(4.37) 

(4.38) 



The coefficients Ci, C2 and C3 can like jFi, JF2 and JF3 be expressed as convergent integrals 
and will be calculated in the next sections. 



5 General Strategy and Technical Tools 

In this section we give a brief overview over the analytical tools involved in the calculation 
of the Feynman-diagrams. These are the measure (subsection 5.1), the procedure to 
extract the residue and the conformal mapping (subsection 5.2). 

5.1 Analytic continuation of the measure 

We first define the explicit form for the integration measure in non-intcgcr dimension D, 
that will be used in the calculations. We use the general formalism of distance geometry 
[37], which has been already used to construct other, but equivalent, representation for 
such measures. 

The general problem is to integrate a function /(xi, ■ ■ ■ , Xn)-, which is invariant under 
Euclidean displacements (and therefore depends only on the N{N— 1)/2 relative distances 
\xi — Xj\ between these points) over the — 1 first points (the last point is fixed, using 
translational invariance) in IR^ for non-integer D. In order to define the integration, let us 
take D > N — 1 and integer. For i < N wc denote hy i/i — Xi — x^ the i'th distance- vector 
and by its a'th component {a — 1, . . . ,D). 
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The integral over yi is simple: Using rotation invariance, we fix yi to have only the 
a — 1 component non-zero. The measure becomes 

d% = Snl dyl{ylf-\ y^ = {yl,0, . . . ,0) (5.1) 
where Sd is the volume of the unit-sphere in IR"^, defined by 

We now fix y2 to have only a = 1 and a = 2 as non-zero components. The integral over 
y2 consists of the integration along the direction fixed by yi and the integration in the 
orthogonal space IR'^^^: 

//•oo roo 
d% = Sn-i / dyl / dyl {yj)''-' , y^ = {yl, yl 0, . . . , 0) (5.3) 
J-oo Jo 

For the j-th point, one proceeds recursively to integrate first over the hyperplane defined 
by 7/1, . . . , and then the orthogonal complement: 

/ roo roo 
d% = n / / dy'j iyT~' ' yj = (yj, ■ ■ ■ , yj, O, . . . , O) (5.4) 

„^ ■ J—oo JO 

The final result for an integral over all configurations of N points is 



I 



n d% = SnSn-1 ■ ■ ■ 11 11 dy« dy] [y])^-^ (5.5) 



j=l j=l \a=l 



This expression for the measure, now written in terms of the N{N — 1) variables y^, can 
be analytically continued to non- integer D. For D < N — 1 this measure is not integrable 
when some y^ ~ 0. For D not integer, the integration is defined through the standard 
finite-part prescription. This means that the measure (5.5) becomes a distribution. 

Let us made this exphcit on the example of = 3 points, following [37, 34]. The 
measure is then 

roo /'+00 roo 

SdSo-i / dyl {ylf-' / dyl / dyl {y^f-' (5.6) 

Jo J-oo Jo 

It is well defined and integrable for D > 1. For D = 1 the integral over y| diverges 
logarithmically at |/| 0, but this singularity is cancelled by the zero of Sd-i and the 
measure becomes 

roo r+co roo r n 

2 / dyl / dyl / dylS{yl) = dy^ dy^ (5.7) 

Jo J-oo Jo JJR JM 

thus it reduces to the measure for two points on a line. For < D < 1 the integral over 
1/1 diverges at |/| 0, but this divergence is treated by a finite part prescription. 

For integrals over N > 3 points, a finite part prescription has already be used for 
D < 2. This will be shown explicitly later. The expression (5.5) is equivalent to the 
measures defined in [37]. 
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5.2 Extraction of the residue 



We now explain how we extract the residue of the pole at £ = for the example of the 
1-loop counterterms. Note from (3.1) that 



2D Jx<L 



2De 



We used the normalization of the measure 



/ = ^ / d^x , Sn = 2 



D/2 



r {D/2) 



(5.8) 



(5.9) 



which was chosen to simplify the calculations (see appendix A). The residue can most 
easily be extracted by applying to (5.8). This yields: 



X 



D-vd 



= L' 

2D 



So the residue of (5.8) is 



We can apply this recipe to the second 1-loop counterterm: 

/ / 

Jx<L Jy 



u ju^j-i Jy<L y 

since it is also proportional to U . We thus have to calculate 



r d 



+ 

.Jx<y=L Jy<x=L 



-d/2 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



We now introduce a general method which is very useful to manipulate and simplify 
such integrals. It relies on (global) conformal transformations in position space and is 
called conformal mapping of sectors. It has first been introduced in [34] , where a geometric 
interpretation can be found. We will explain the method on a concrete example and then 
state the general result. 

Let us consider the second integral on the r.h.s. of (5.13) 



/ ( 

Jy<x=L ^ 



2u I 2u 

X + y 



-d/2 



L / — X / — y [x + 
Jo X JO y ^ 



,2u 



-d/2 



Now two changes of variables are performed: The first one 

X ^ X , X — xyL~^ 

leads to 



6{x - L)Q{y < x) 

(5.14) 

(5.15) 



Jo X Jo y 



{f + L''') '^'^ 5{xyL-^ - L)Q{L < i) (5.16) 
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The second one 
finally gives 



(5.17) 



Jo X Jo y ^ ' 
Replacing the second integral in the r.h.s. of (5.13) by (5.18) gives 

^1; ( I — - f f 7^" + »^") — »>-'*(^-^> 

(6.19) 

Now one distance (here y) is fixed, whereas the integral over the other distance (here x) 
runs from to oo. The former constraint max(a;, y) = L has been transformed into the 
factor max(a;, |/)"^ times the constraint y — L. 

Before generalizing this formula, we shall show how it can be used in practice. The 
residue in - (which determines the corresponding 1-loop counterterm) is given by the 
simple formula {d^{D) = 4L)/(2 - £>)): 



°° da; 



X 



(f'' + l) 



-dc{D)/2 



D 

2-D 



2-D 



r(^) 



(5.20) 



The subleading term can analogously be calculated by expanding (x^'^ + 1) dc{D))/2 
max(x, yY in e. We obtain the convergent integral representation 



L 



-X 



D 



X 



2-D 



\n.{x ^ + 1) — ln(max(a;, 1)) 



(5.21) 



This method extends to the integrals which appear in the counterterms associated to 
the contraction of any number of points. In general we have to compute integrals over 
N{N — 1) distances x,y,- ■ •, of the form 

I{s) = f /(a:,y,...) (5.22) 

Jinax(x,y,...)<L 

with an homogeneous function / such that the integral has a conformal weight (dimension 
in L) K. For the integrals which appear in n-loop diagrams, this weight is simply 

K — ne (5.23) 

The integral over the distances is defined by the D-dimensional measure (5.5). The residue 
is extracted from the dimensionless integral 

J{e) = L-'^L^Iie) = neL-^I{e) ^ L [ /(x, y, . . .) max(x, y, . . .)-'^ 

OL/ Jm.ax{x,y,...)=L 

(5.24) 
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The domain of integration can be decomposed into "sectors" , for instance 

{■■■<y<x = L}, {■■■<x<y = L} , 



(5.25) 



and we can map these different sectors onto each other by global conformal transforma- 
tions. For instance we can rewrite the integral (5.24) 

J{e) = L f{x, y,...) max(x, y, . . .)" = L / f{x, y,...) max(x, y, . . .f (5.26) 

Jx=L Jy=L 

The constraint on the maximum of the distances is replaced by the constraint on a (arbi- 
trarily chosen) distance. 

This mapping of sectors is one of the basic tools used in the following to explicitly 
calculate the 2-loop-diagrams. 

6 Coupling constant renormalization, first graph 
6.1 The count erterm 



Figure 6.1: The distances in (6.2) 



We are now going to calculate the first diagram of section 4. It contributes to the 
coupling-constant renormalization in 2-loop order and is 



3 



(6.1) 



With the distances labeled as in figure 6.1, the first MOPE-coefficient is: 



• • 



^ (^\/^2^^Td2^+ VP^^T^+ V^^r^T^^) [yh'^^ + e'^^ + ^Jc^^ + - Va^'^ + d?""^ 



(6.2) 

The derivation of this expression can be found in appendix B. Recall that the factorization 
of the divergences, which are subtracted in (6.1), is: 



(6.3) 
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We therefore represent the MOPE-coefRcients associated to the second term in (6.1) as: 





-d/2 



-d/2 



(6.4) 



There are 3 sub divergencies subtracted due to the 3 possible contractions and a factor 
1/2 due to symmetry, i.e. due to the nested contraction. 

We can now proceed to calculate the diagram (6.1), which is of order 1/s and not of 
order 1/s^ as the single terms. To do so, let us consider the integral I{L), which is defined 
as the integral of the MOPE coefficients, with all mutual distances appearing in figure 6.1 
restricted to be smaller than L: 



I{L) 



a,b,c,d,e,f<iL 






(6.5) 

Note that /(L) is no^ exactly equal to (6.1). However we shall show below that the leading 
term, i.e. the pole in 1/s in which we are interested, is the same. 
I{L) has the following Laurent-expansion: 



I{L) = + 0{s^: 



(6.6) 



We now apply the operator to I{L) to extract the residue a in 1/e. We obtain 
an integral similar to (6.5) with the constraint that max(a, b, c, d, e, /) = L. Using the 
trick of conformal mapping explained before, we can rewrite this as an integral with the 
constraint that one of the distances (for instance a) is equal to L and that the other can 
vary freely. We thus obtain the integral 



J(L) 



' lllllb cdef ^' ^' ^' ^' " 



■2e 



X 



X 





(6.7) 



We know that the only divergencies, which lead to poles in 1/e, could appear in (6.7) when 
a pair of distances, i.e. (a, d), {b, e) or (c, /) simultaneously tends to 0. But by construction 
these divergences cancel between the first and the second terms of the integrand. So, the 
integrand in (6.7) has only integrable singularities at short distances and we can perform 
the limit d dc, i.e. £ — > inside the integral in order to compute 



(6.8) 



This convergent integral can be calculated numerically by the methods developed in [34] . 
We shall explicit this calculation later. 
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The next problem we have to treat is that in (6.5) we subtracted a counterterm that 
we had chosen by convenience. It is however not the counterterm in (6.1), which is fixed 
by our renormahzation prescription at 1-loop order. The latter can be written as 

The difference between (6.9) and the counterterm in (6.5) is that the integration over c 
and / is not restricted. To evaluate the residue for the true counterterm, we apply the 
operator to (6.9) and get an integral with the constraint that max(a, 6, d, e) = L. We 
now use the trick of conformal mapping to transform this integral into an integral with 
the constraint a = L as in (6.7). This yields: 

^IIIIL d ^{^''^ + d"'y'^\b"' + e"'y'^'im.x{a,b,d,e)-'' (6.10) 

This expression has to be compared to the corresponding counterterm in (6.7). The 
difference between these two counterterms is 

///// (^''^ + '^'^) (^'' + ^'') H""^"' - •^)"'') 

J J J J J h^c^d^G-^f 

(6.11) 

A priori this term is 0{e^) and might contribute to the residue. We state that this term 
is 0{e), thus subdominant and does not contribute to the counterterm. To prove this we 

set L = 1 and wc develop max(. . .)^^^ in powers of e. This expansion is licit over the 
domain of integration, since max(. . .) > 1 and since the integral is convergent at large 
distances. This implies that (6.11) is equivalent to 

2e jjjjj (0"" + d''^) (b'"' + e^") '^''^ [ln(max(a, 6, c, d, e, /)) - ln(max(a, 6, d, e))] 

(6.12) 

when £ — > 0. Poles of the integral in 1/e might appear if 6, e — > simultaneously. In that 

case however c ^ a and d f so that the difference in the last factor in (6.12) vanishes. 
The last integral in (6.12) is therefore convergent, implying that (6.11) = 0{e) and that 
/(L) + 3/2C1 = 0{s^). Therefore the use of (6.5) instead of (6.1) to compute the residue 
of the pole is justified. 

Such a phenomenon is not peculiar to this diagram. In the other diagrams that we 
shall calculate similar simplifications occur when dealing with the counterterms, which 
allow to take the same constraints over the distances for the 2-loop diagram and for the 
counterterms associated with 1-loop sub divergences. We call this property the "2-loop 
miracle". (A generalization to the calculation of higher order diagrams is possible.) 

6.2 Numerical calculation 

We now want to calculate the residue J(L), given by the integral (6.7), numerically. This 
calculation will be performed for values of D in the interval 1 < D < 2. As we already 
discussed, the integral representation (6.7) suffers from additional divergences for D < 1, 
which come from the fact that the measure over three points becomes a distribution, and 
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a 



Figure 6.2; The half-sector used for the numerical integration 



which must be treated by a finite part prescription. This will not be done here, since it 
turns out to be sufficient to calculate the diagrams for 1 < D < 2 in order to have good 

estimates for the critical exponent of membranes. 

The calculation is considerably simplified by using the symmetries of the diagram 
and by integrating over one of the equivalent sectors only. First of all, one can suppose 
that max(a, &, c) > max(d, e,/). Furthermore an ordering of a, b and c is introduced: 
a > b > c. This gives rise to a symmetry- factor 2 x 6 = 12. In order to eliminate all 
divergencies, also d, e and / have to be ordered. We assume that d > e > f but then 
have to sum over all permutations of d, e and /. Denoting by 



g{a, b, c, d, e, /) = 



as in (6.2) and (6.4), we get 





(6.13) 



J{L) 



6=0 



a=l,b,c,d,e,f 



g{a, b, c, d,e,f) 



^2 ja=i>b>c ^' ^' ^' /) + 9{a, 6, c, e, /, d) + g{a, 6, c, /, d, e) 

a>d>e>f 

+9{a: b, c, /, e, d) + g{a, 6, c, e, d, f) + g{a, 6, c, d, /, e) (6.14) 



We can furthermore divide the remaining integral into 2 sectors where either 6 < e or 
b > e. There will always be an integrable divergence in the smaller distance which has 
to be treated by an appropriate variable transformation. As in [34] we parametrize the 
integral over b as: 

Sd-i 



Sj^ J — oo 

1, b 



dbi 



db2 b2 ^ 



(6.15) 



bl 



restricted to the domain where b < c < a: There are singularities for 6i ^ and for 
62 0. They are disentangled by switching to radial and angular coordinates. In these 
coordinates the singularities can be eliminated by the following parametrization: 

1_7tSd^i '-^ 

lb ^ ^ - 



D-aD-12 S 



D 



df3 pT^^ sm{af-^ x 



X 



C dtb^Qil-b^Qic-b) , 
^0 



(6.16) 
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where 




(6.17) 

(6.18) 
(6.19) 

^62 + i-26cos(a) (6.20) 

The change of variables from the angle a to /?, equation (6.17), generates a factor 
f]{'^-D)/{D-i) -^y^yq}^ exactly cancels the singularity for a — > from the measure, i.e. 
sin(Q;)^~^. The other variable-transformations are constructed in a similar way. 
Furthermore there is the integral over d and e. 



Idle 



I with e' = ^ (6.21) 

d Je' d 



Therefore we write 



/ dDd=l-['dr d = r'/^'' (6.22) 
Jd<a 2D Jo 

The remaining integral over e' is parametrized in analogy to the integral over b as: 

f dt{e'yQ{d-f)Q{f-e) , (6.23) 

where 



7 =-5^5-1 (6.24) 
2 



e^dsD-T (6.25) 

_ f , e < 6 
^~l2D-2, e>h 



(6.26) 



/ = ^d2 + e2-2edcos(7) (6.27) 

With these variable transformations, the integrand is bounded. This does not mean that 
the numerical integration is easy. The main problem is that the integral is localized in 
some small domain, in which the integrand (with all the factors of the measure and from 
the variable transformations) is about 100 or 1000 times its mean-value, whereas it is much 
smaller in large domains of integration. A genuine adaptive Monte-Carlo (AMC) routine 
has to be used. It is described in our earlier publication [34]. The idea is to divide the 
domain of integration into subboxes and to try to integrate each subbox using standard 
Monte-Carlo (MC) integration with few sample-points. The MC-routine gives an estimate 
for the integral and for the standard deviation from this value. If the latter is too large, 
the box is divided into smaller subboxes and the procedure is repeated. This algorithm 
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•Q-B □ m □ a □ m □! 



1.6 1.7 1.8 1.9 2.0 



£=0 




D 




1.00 


1.5 


1.01 


1.441 ±5 X 10^3 ±4 X 10-3 


1.02 


1.399 ± 5 X 10-3 


1.05 


1.268 ±5 X 10-3 


1.10 


1.065 ±5 X 10-3 


1.15 


0.876 ±4 X 10-3 


1.20 


0.702 ±4 X 10-3 


1.25 


0.543 ± 3 X 10-3 


1.30 


0.399 ± 2 X 10-3 


1.35 


0.275 ± 2 X 10-3 


1.40 


0.176 ±2 X 10-3 


1.45 


0.1022±7 X 10-^ 


1.50 


5.16 X 10-2 ±4 X 10-^ 



D 




1.55 


2.16 X 10-2 ±2 X 10-3 


1.60 


7.07 X 10-3 ±6 X 10-^ 


1.65 


1.61 X 10-3 ± 1 X 10-^ 


1.70 


2.17 X 10-^ ± 2 X 10-^ 


1.75 


1.23 X 10-^ ±2 X 10-^ 


1.80 


1.55 X 10-^ ±2 X 10-9 


1.85 


1.07 X 10-1° ±2 X 10-12 


1.90 


4.72 X 10-" ±6 X 10-19 


1.95 


3.99 X 10-36 ± 6 X 10-36 


1.975 


2.81 X 10-^* ± 3 X 10-^6 


1.98 


2.33 X 10-93 ± 2 X 10-95 


1.99 


6.89 X 10-189 ± 1 X 10-189 


2.00 






Figure 6.3: Numerical results for the diagram (6.7). The first error in the table is the statistical error, 
the second the systematic error. The latter is only given if it can not be neglected. 



was implemented recursively using the computer- language C. For details cf. [34]. The 
AMC-routine gives an estimate of the integral and of its statistical error. A systematic 
error also appears, which is more difficult to estimate. It comes from the domain of small 
angular variables, where the numerical precision of the workstation is no longer sufficient 
and it appears to be the most important for — * 1, as can be seen from (6.16). It can 
be estimated by counting exceptions of the floating-point unit and comparing it to known 
integrals. For (6.7), the systematic error is negligible. Its treatment will be discussed for 
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the diagrams, where its contribution is important. The results for (6.7) where obtained 
within some hours on a workstation and are hsted in figure 6.3. The numerical results 
nicely fit with the analytical value for D — 1, discussed in the next section. 

We checked the numerical integration routine. This check is provided by integrating 
the function 

g{a, b, c, d, e, f) = {a^ + b'^ + + (f + + fY^^ (6.28) 
which has conformal weight k = and by comparing to the analytical solution 

/ gia, b, c, d, e, /) = . (6.29) 

Ja=i,b,c,d,e,f ' ' ' ' '•'^ 8-3^r(2D) ^ ^ 

This test ensures that the integration can be done for 1 < D < 2. 
6.3 Analytical calculation for D — ^ 1 



— •- 



c 



Figure 6.4: Connection between the MOPE-coefScient and the diagrams in polymer theory 

To check the numerical calculation, we want to evaluate (6.7) for D — 1. This calcu- 
lation can be performed analytically. 

Some subtleties have to be taken care of in order to understand the calculation which 
follows. We remarked that in (6.7) the factor max(a, 6, c, e, Z)"^*^ could be dropped 
without changing the result as the contribution to the integral is finite in any subdomain. 
For the following calculation we shall drop this factor but shall not take the limit d — > d, 
as we only can calculate the diagram and its counterterms separately. Single terms will 
thus have divergencies in 1/e, which arc treated by a finite part integration prescriptions, 
and which have to cancel at the end. Through this change only the sum of all terms but 
not each single term has a meaning in the limit £ — > 0, that we take at the end. 

For D — > 1 the measure locahzes on a fine. This introduces different orderings of 
the distances, which are topologically inequivalent. Two types of diagrams appear for 
D ^ 1. There arc cither the "untwisted" diagrams (figure 6.5) or the "twisted" diagrams 
(figure 6.6). We remark that the second diagram in figure 6.5 is untwisted as it can 



Figure 6.5: The 2 untwisted diagrams 

be transformed into the first one by simply exchanging the orientation of the lower line. 
With the same reasoning one deduces that the 4 diagrams in figure 6.6 are topologically 
equivalent. 
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Figure 6.6: The 4 twisted diagrams 



Let us start to calculate the untwisted diagram (without the counterterm) . We note 
that d = A-2£. We find 



b=l 


a 




d 



-* » *- 



= / / [ {ab + ae + bd + ed)-"^/^ 

J a J d Je 

= / / /(a + ci)-'^/2(6 + e)-<^/2^0 

J a J d J e 



(6.30) 



The fact that the last integral is zero is a well known property of the finite part integration 
of a homogeneous function (here j^{a + d)~'^/'^—Qi). 
The counterterms are: 



a J d J e 



{a + dy'^b + e)-^/^ = Q (6.31) 
/ / [(a + dy'^/^c + f)-'/'^ [ [ f{a + d)-'/'il + a + d + e)-'/' 

JaJdJe JaJdJe 



11 1 X 



a J d J e 



{b + e)"^/^ (c + /)-<^/2 
So together these terms add up to 



a Jd Je 



^l + e)-''/'(l + a + d + e)-'^/' 



6 



1 \ 3 

- ^ counter terms j ~ ^ + ^(^) 



(6.32) 



(6.33) 



(6.34) 



The factor 6 is the combinatorical factor. 
The untwisted diagram is given by: 



b=l \ a ; 



f <> 



= / / / (ab + af + ad + bd + df)"^/^ 

J a JdJf 

^111 {a + af + ad + d + df)-'^/^^]-- + l + 0{e) (6.35) 

JaJdJf 8£ 8 



The counterterms arc: 

(a + d)-'^/^ (i, + e)-'^/2 



a Jd Jf 



a Jd Jf 



{a + d)-'^/^l + d + f)-'^/^ 



11 1 



(6.36) 
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{a + d)-''/^c + fy^^ = / / / (a + d)-'^/2(l + a + /)-'^/2 

laJdJf JaJdJf 

^ 0{e) (6.37) 



a Jd Jf 

= l + 0(e) 

These terms add up to 

12 



8£ 4 



/ / / (l + / + d)-<^/2(i + a + /)-'^/2 

J a Jd J f 



1 \ 3 

— - counter terms j = - + 0{e) 



(6.38) 



(6.39) 



The factor 12 again is the combinatorical factor. 

The final result which has to be compared to the numerics thus is (cf. (6.7)): 



J{L = 1) 



_ 3 

D=l ^ 



(6.40) 



7 Coupling constant renormalization, second graph 




Figure 7.1: The distances in (7.1). 



7.1 Derivation of an analytical expression 

The next diagram which has to be calculated is (4.37): 



2\'^ 




(7.1) 

The second term on the r.h.s. subtracts the subdivergence, when first the single dipole on 
the r.h.s. of ^^^^^ tf^'') is contracted. In this case the MOPE is: 



Note that the MOPE (7.2) has a tensorial structure and that a counterterm like 



(7.2) 



(7.3) 
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would subtract the pole-term from the integral but would not be sufficient to make the 
integral convergent. With the distances as noted in figure 7.1 the MOPE-coefficients are: 



2. .^2u^^2u 



d 



2(,2v I 2i/\-(i/2-l r-L'{d+2) 



/-C") {bfb-" - cjc-") 



-d/2 

(7.4) 

2 



+ ^2.yd/2-lj-.(d+2) (^^j^-D _ dfd-^y 

(7.5) 



where two equivalent formulations were given. 

For scalar-products we always use the convention that the vectors involved start from 
the same point, i.e. we define (cf. figure 7.1) 



bf-^{b' + f-d') 



(7.6) 



A relevant counter-term, which we did not explicitly write in C2, appears too. It is canceled 

by ^ ^ 

' ^ ' ' ^ (7.7) 



We state and will show below that 
I{L) 



a,b,c,d,eJ,g<L 



^ )e(/ < max(a,d,e,t/)) 



(7. 



where we used the second version of the counter-term in (7.5). 

To compute the residue of the single pole, we apply Ld/dL on this expression and 
map onto c — L = l: 



(-^) ^ -^^-^(-^) ^ ////// max(a, 6, c, d, e, /, g) 

OLi J J J J J J c=l\a,bA,eJ M 



; a,b,d,e,f ,g 



-d/2 



-2e 



(y2^ + a^'^y'^/^f-'^'^ 



-t'ig'^ + a'T'^'-'n^'^'^ (efe-^ - dfd-^Y Q{f < max(a, d, e, g))^ (7.9) 

We check that this expression is integrable everywhere, as is indeed the case and can be 
seen by a (generalized) Taylor-expansion. 
We then have to explain that 



/ 

Ja,d,e,g<L;f \ 



^ ) e(/ < max(a,d,e,^)) 
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I e(/ < max{a,d,e,g)) 



(7.10) 



The first equafity is due to symmetry. The second stems from the so-called nested inte- 
gration. We see this exphcitly as follows: 




) e(/ < max(a, d,e,g)) 
I max(a, d, e,^g)\ 



l,e,g<L \ 

where we used the fact that for any / 



(7.11) 



(7.12) 



The whole integral scales like L^^, so that we can apply to (7.11) without changing 
it. Doing so we obtain that (7.11) equals 



L 

2,£ Jniax.(a,d,e,g)=L 



I max(a, e, g)' 



X 



^) (7.13) 



The factor ^ '"^(M.e.g) ^ jg equal to unity and the integral is our usual expression for the 



residue of 



; . Taking all this together we obtain 




(7.14) 



This proofs the disered result. The reader can verify that this is a general feature of these 
so-called nested contractions. If the largest distance in the subdiagram is restricted to be 
smaller than the largest distance which rests after complete contraction of the subdiagram, 
then the so restricted integral is 1/2 times the product of the subdiagram and the diagram 
which rests after contraction. This demands of course that both of them scale like U , 
where L is the IR-cutoff. The reader will be able to generalize this rule for a scaling with 
other exponents, which would be necessary in higher order calculations. 

The subtle point which we still have to check is that the changes in the domain of 
integration of the marginal counterterm from (7.1) to (7.8) do not change the residue, i.e. 
that again the "2-Ioop- miracle" appears. Analogously to section 6.1 we write down the 
difference, apply Ld/dL on this expression and map onto c = 1. We obtain 



d 



Illllla bdefg {^^^^^' ^' ^' ^' ^' ~ niax(a, d, e, /, g) 
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We would like to develop (max(a, b, c, d, e, /, g)^"^^ — max(a, b, c, /, g)^^"^) for e small and 
show that the corrections are of order e. This might be wrong, if and only if this expression 
does not vanish at points, where the integral has a pole. For / ^ it vanishes. The 
limit a, (7 — >■ is a bit more subtle as the difference does not vanish. However in this case 

{efe~^ — dfd~^^ is of order so that no pole in the integration over a and g appears. 

Let us now perform in (7.9) the limit d dc, i.e. e ^ 0. Then we would like to 
integrate over g analytically. This is not possible due to the ©-function. We therefore 
modify this constraint from 0(/ < ma.x{a,d,e, g)) to 0(/ < d). Note that the modified 
counterterm still successfully subtracts the marginal subdivergence. We obtain: 



"^^ ^ lllllla,b,d,e,f,! 



-dc/2 



-^v\g^^ + a2-)-'^=/2-i/-(rfc+2) ^^gj^-D _ dfd-DyQ^f < (7.16) 

Of course this change affects the result and we shall calculate the difference in section 
7.7. The modified counterterm has another useful property, which justifies its choice: The 
©-function is not affected by the R-operation discussed in the next section. 
Now the integration over g can be performed. The result is: 

J(L) =///// F{a,b,c,d,eJ) (7.17) 

J J J J Ja,b,d,e,f 

with 

, 2 



1 r 

F{a, b, c, d, e, /) 



D 

2-D 



2-D 



2\ 



-y^a-V-'-^e(/ < d) (efe-^ - d/d-^f } (7.18) 
7.2 Improvement of the measure 

For c—1 the measure, given by (5.5) and (A.l), simplifies to an integral over the vectors 
a and /: 

% / / da^at' d/i / d/^ / df, f^-'F{a,b,c,d,e, f) (7.19) 

0£) J —oo Jo J —oo J —oo Jo 

For D < 2 the measure defined in (7.19) is a distribution and suffers from a relevant 
divergence for /3 — > 0. Geometrically these are configurations, where the tetrahedron 
spanned by a, 6, . . . , / has volume 0, i.e. is restricted to a plane. A finite part prescription 
has to be applied in order to make the measure finite. This was first discussed in [34] . 
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One may think of implementing this prescription by subtracting the singularity. This 
method however imposes at least numerical difficulties. It is better to eliminate the singu- 
larity by a partial integration with respect to fs, which is mathematically equivalent [34]. 
As only d, e and / depend on /s, the integral 

POO 

f.p./ d/3/3^-V(a,6,c,rf,e,/) (7.20) 



can be converted to 

1 roo f) 

= 2^/0 d/3/3''"'Ri^(a,6,c,d,e,/) (7.21) 

where R is defined via 

The strength of the divergences for a — > or / — > is unchanged. It is important to 
remark that this trick cannot be used to eliminate the relevant divergences when a — > or 
/ — s> 0. It works for the divergence in /s, because the integrand does not directly depend 
on /a but on rf, e and /, which themselves depend on /a. So the derivation of F with 
respect to /s does not produce a factor I//3 but factors l/d, 1/e and 1//, which are not 
singular for /a — > 0. Explicitly: 

1 

RF(a, h, c, d, e, f) = — — — -A- x 

2- D r (y^) 



D 



D f-2D-2 



+2Da-^/ 

+(2 + D)^a-'f-'-^e{d > /) {efe-^ - drd-^'f 
-Dua-^f-^-^e{d > f) [efe-^ - dfd'^) x 

X (2e-^ - 2d-^ - Defe-''-^ + Ddfd-''-^) | 



(7.23) 



Note that R0(/ < d) = whereas R0(/ < max(a, d, e)) gives contributions proportional 
to e.g. 6{a — f) which had to be treated separately. This justifies our choice of the modified 
bound of the counterterm in (7.16). 
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7.3 Parametrization of the measure 



The main singularities for small distances appear for a or f small. We therefore want 
to parametrize the measure with the help of these distances. The divergences for small 
volume of the tetrahedron spanned by a, . . . , /, a2 ^ or ^ shall be treated by a 
parametrization in angles as by this way small distance and small volume singularities 
are best disentangled. We have chosen the parametrization indicated in figure 7.2. One 



c=l 

Figure 7.2: Parametrization of the tetrahedron 

triangle is spanned by c and a with an angle /? between, another by c and /, where the 
corresponding angle is a. The angle between the planes spanned by these two triangles 
is r. The distances as functions of a, / and (3, a, t are: 



6 = Y + 1 — 2a cos /? 

e = ^/2 + 1-2/ cos (7 (7.24) 

d — (a cos (3 — 1 + f cos (t)^ + (a sin /3 — / sin a cos t)^ + (/ sin r sin 

The integrals over a and / run from to oo, the integrals over /3, a and r over the interval 
[0, tt] . As wc do not want to map all the points which are far away, we have to find a 
reparametrization of the measure which behaves for a ^ like a'^ and for a — oo like a*^, 
by this way eliminating the principle divergences. If u is equally distributed we can use 

1 1 

a^uD-'r{l-u)D-u> . (7.25) 
The integral over /? will be parametrized as 

71 



^ ! _j_ ~ (7.26) 

t: - -{2-2a)D-i a > 0.5 

The integral over / (note that the factor /| came from the partial integration with 
respect to /a) 

JdhdhdfsftVi) (7.27) 
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can be written as 

/ d/r-i/ da{smafj drismrf-'if) . (7.28) 
We change variables from f to v: 



f = vD-'r{l-v)D-u^ (7.29) 
Furthermore we choose in the same spirit as for (3 

TT 



a 



^{2r])D+i ?7<0.5 



TT 



TT- -(2-277)^+1 77 > 0.5 



(7.30) 



and 



|(2C)^ C<0.5 

7r-|(2-2C):^ C>0.5 

So the complete integral over four points is: 

Sd-iSd-2 

si {2-D){D-l)D{D + l) 

'■1 n / 1 1 1 1 \ . _ „ „ 



(7.31) 



/ dua^ I + 1/ da min(2Q;, 2 - 2q;) i^-i (sin /3)^~^ 

Jo \D- 7 u uj- Dl- uj Jo ^ ' ) \ f^) 

C dvf (—^1 + ^--^) /'dr7min(2r7,2-2r7)^(sin(7)^ 

f dC min(2C, 2 - 2C)^(sinT)^-V'Ri^(a, 1, c?, e, /) (7.32) 
Jo 

Another way of parametrizing consists in replacing the integral over the vectors a and 
/ by the integral over the vectors a and d. This parametrization is especially useful to 
eliminate divergences, when a and / simultaneously go to infinity. It will be used for the 
integration over one of the sectors in the next section. The new formulas are given here, 
a prime indicating new angles and distances as can be deduced from the figure, a' and 
t' obey the same relations as a and r. For d' we use the same variable transformation as 
before for /: 



d' ^vD-i{l-v)D-u, (7.33) 

The other new distances are: 

e' = Va^ + d'^-2ad'cosa' (7.34) 

/' = ^J {d' sin a' sin t'Y + [d' cos a' — a + cos /3)^ + (sin (3 — d' sin a' cos t'Y (7.35) 

For the integrand (7.23), the exponents 7 and uj are found by performing a (generalized) 
Taylor expansion: 

7 = (7.36) 

uj^2D (7.37) 
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c=l 

Figure 7.3: Alternative parametrization of the tetrahedron 

7.4 Decomposition into sectors 

Although the measure absorbs the principal singularities it cannot handle all of them. 
There remains e.g. a singularity for 6 — > and e — > 0. Two methods may be applied 
to handle the remaining integrable singularities. The first consists in using the second 
measure of section 7.3. The second is to map again some parts of the domain of integration. 
Thereby, we face the problem that the measure is no longer symmetric in the distances, 
as we have changed it in order to eliminate the relevant singularity for /a 0. In order 
to restore this symmetry, we rewrite the integral (7.32) as 

Sd-iSd-2 tt^ 



SI (2- 

dua'^^^ { — - + — ^ ^ ) / da min(2Q;, 2 - 2Q;)^(sin/3)^ 



I 
I 




2-D 



with 



dv ( — ^ — - + — ^ ^ ^ ) / dr] min(277, 2 - 2r]) d+i (sin a)^ 
yD — ^ V ■ -.1/ 

dCmin(2C,2-2C)"^(sinT)^-^r(a,6,c,(i,e,/) , (7.38) 



T{a, b, c, d, e, /) = ————RF{a, b, c, d, e, /) , (7.39) 

(2A(a, b, c)y 



and where A(a, b, c) = |ac sin(/5) = {a + b + c){a + b — c^{b + c — a){c + a — b) is the 
area of the triangle spanned by a, b and c. This is, except for the geometric prefactor, 
the invariant measure in D + 2 dimensions. The integrand now is conformal invariant, as 
follows directly from equations (5.22) to (5.26). 
The sectors are decomposed as follows: 

1. (a < 2 or / < 2) and {b > I or e > I) 

This sector is convergent: F{a, b, c, d, e, /) is integrated directly, using the simple 
measure (7.32). 
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Figure 7.4: Numerical results for (7.23). The first error is the statistical error, the second an estimate 
for the correction due to the systematic error, which becomes important for Z) — > 1, cf. the text. In the 
plot, the boxes are the uncorrected, the crosses the corrected results. 



2. a > 2 and / > 2 

The measure (7.32) does not eliminate the singularity, when both a and / simul- 
taneously go to oo. The easiest way to integrate this sector is to use the second 
measure (7.33) ff. of section 7.3. 

The divergences of the integrand could also be eliminated by a mapping. This 
however induces new singularities due to the measure (the term l/(2A(a, 6, c))2 in 
T, equation (7.39)). This would not be the case, if we had not been forced to use 
the trick of integrating by parts the measure. 

3. 6 < I and e < | 

In this sector the mapping can be used successfully: a has to be exchanged with h 
and e with /. We get T{b,a,c,d,f,e) with a < | and / < |. This is integrated 
using the measure (7.38). 
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7.5 Numerical calculations 



The numerical calculations are difficult. We refer the interested reader to the discussion 
in section 6.2 and for more details to [34], section 6. Here we only give the result, see 
figure 7.4. The extrapolation for D — > 1 is consistent with the analytic result 2, found in 
subsection 7.6. 

We however want to discuss our estimate of the systematic error. It is is a general 
phenomenon that our algorithm fails to correctly integrate for integrals over 4 points in 
the limit D — > 1. This is due to problems in the domain 02 — > and /2 ^ and thus 
occurs even for integrals which are well-behaved at small and large distances. An example 
of such a function is 



F(a, 6, c, d, e, /) = (a^ + + + + + /2)-30/2 
which is analytically integrated to give 



1 (-Y 

shr(^) V4 



(7.40) 



(7.41) 



The ratio of the numerically and analytically calculated values is displayed in figure 7.5. 
As these factors should be independent of the integral which has to be performed, we use 
them to correct the numerical results. 
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Figure 7.5: Numerical results for the integral of (7.40) 



7.6 The limit D ^1 

In the limit D — > 1, (7.17) can again be calculated analytically. This calculation is 
interesting as it reveals the connection to standard polymer theory and the fact that 




* j decomposes into three topologically different and non-equiva- 



lent diagrams. As in equation (7.17) we keep c = 1 fixed. By a direct calculation it can be 
verified that the measure indeed reduces to an integral over a line. On this line two points, 
the endpoints of c, are already fixed. Then there are 12 different possibilities to distribute 
the last two points. They still can be separated into four topological inequivalent classes A, 
B, C and D, cf. figure 7.6. These are the four standard diagrams arising in polymer theory. 
In each of these classes the line with c = 1 may be chosen to be the line connecting (12), 
(14), (23) or (34). Readers more familiar with Feynman diagrams arising in the framework 
of a scalar field theory may recover the three corresponding diagrams after a de Gennes 
transformation [40] . They contribute to the renormalization of the (/j^- interaction aX d — A 
and are represented on the r.h.s.. Diagrams in one class can be mapped onto each other 
by the now well-known mapping of sectors. One subtlety however has to be taken into 
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A = 



B = 



C = 



D 



2 4 



• • 



1 1 3 




Figure 7.6: The 4 topological inequivalent classes A, B, C and D and the equivalent diagrams in the 
n ^ limit of scalar (^'* theory. 



account. The marginal counterterm in equation (7.18) is not invariant by this mapping. 
In order to perform the integration only over one sector in every class, the symmetrized 
version will be used: 

F(a, 6, c, d, e, /) = r\af -\{h-c + e- df)'' - a" V"' 



The diagrams give: 



16 



(e(/ < c) + e(/ < b)) (bfb-' - cfc'] 

+ (e(/ <d) + e{f < e)) {efe-' - dfd-'f 



POO POO 



1 r°° 
4 



1 11 



/(a/-l) aP 16 aV 



e(a-l-/) 



- - — ln(2) 
4 16 ^ ^ 



/ da / d/— - — -0 



1 /""^ 
4 



1 Z'^*-' 

4 



/■OO POO 

dd df 

Jo Jo 



1 



ni+d) {i+f+d)f 



(7.42) 



(7.43) 



(7.44) 
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4(l + / + d)2/ 



1 



(2 + e(/<i) + e(/<d)) 



(7.45) 



1 r<x> ff-i 1 



1 

4 



(7.46) 



Taking care of the combinatorical factor 4 for the sectors A and B and of the factor 2 for 
C and D, the final result is: 



+ 4 



+ 2 



+ 2 



(7.47) 



7.7 The correction for the unusual marginal count erterm 



We recall that J{L), given in equation (7.16) and which we calculated numerically, was not 
exactly the counterterm J{L) but was modified in order to simplify the calculations. We 
still have to calculate the difference J{L)—J{L) , which also contributes to the counterterm: 

J(L) - J(L) = ^ / (a- + g'Y"-' r^'^'^ {hfh-^ - cfc-^f x 

Z J c—l, a,b,d,e,f,g 

X [e(/ < c) - e(/ < max(a, b, c, g))] max(a, b, c, f, g)'^' (7.48) 
First of all the integral over / is performed. Since 

/^^^ {bfb-'' - cfc-^y 

we get 



D V 



(7.49) 



J(L) - J(L) = ^ / (a- + ,-)-^/^-^ (6V- - CC-) X 



2-D 



X - [c^ — max(a, 6, c, gY] max(a, b, c, f, g) 



-2£ 



/ (a'^ + g'^y''"'^'Cbb-^-cc X 

■/a=l; b,c,g V / V 

X [ln(c) - ln(max(a, 6, c, g'))] + 0{e) 



(7.50) 



as: 



bmce — cc is symmetric under the exchange of b and c, this can still be written 

2 — Df / o,, n,.\-dc/2-l 

ln(6c) - ln(max(a,6,c,5()) + 0(e) (7.51) 



J(L) - J(L) = ^-^ / (a^^^ + ^2-) "'^'^ ' [h^-^D ^ ^2-2D ^ _ ^2 _ c2)5-^c-^)^ 

2 Ja=l;b,c,q V / V / 



X 



■1 



X 
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Figure 7.7: Numerical results for equation (7.51). The error is ±10 relative. 
Numerical integration 

The domain of integration lias to be split into the two sectors where either 6 or c is the 
smallest distance. These two sectors are equivalent, so the integral will be performed over 
b < c only. The integral is furthermore split into the integral over the radial and the 
angular coordinate. We use the following change of variables: 



b = 1^^-76(1 - 



"6 



-D < 2D 
2 



lb ^ 2D -2 



/5 = 



TT ^— 

-{2a)D-i 



t: - -{2-2a)D-i 
= ^62 + 1 _2fecos(/3) 



a < 0.5 
a > 0.5 



(7.52) 
(7.53) 
(7.54) 

(7.55) 
(7.56) 
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This gives 
~~ D-1 S 



D JO 



da min(2Q;, 2-2q;)o-i sin(/3) 



£1-2 



du 



The integral over g is independently parametrized as: 



\D — ^}jU ujb — D 1 — u 

(7.57) 



g — tli>-79 (^1 _ D-u,g 

7p = 



a;g = l + L><2 + L> 



This imphes 



dt; 



11 11 



-- + 



D 



ig Jo \D — ^gV Ug — D 1 — V J 

The results of the numerical calculations are given in figure 7.7. 
The limit D 

For D —1 the integral can again be performed analytically. We get: 
dg / d6(l + ^)-3(ln(6)-ln(max(l,5))) 



(7.58) 
(7.59) 
(7.60) 



(7.61) 



= / dg (1 + g)-' db Hb) - dgil + g)-'Hg) 
Jo Jo Ji 

= -\ - ^ln(2) = -0.5965735903 
4 2 ^ ^ 



(7.62) 



8 Complementary Contribution for the Renormal- 
ization of the CoupHng Constant 



We still have to calculate C3, equation (4.38): 



U) 



The MOPE-coefficient is: 
The non-trivial diagram in (8.1) is: 



(8.2) 



J a,b,c,g<L y ~- 



(8.3) 
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Apply Ld/dL and map onto a = L = 1: 





max(a, 6, c, g)' 



X max(a, 6, c, ^f)^*^ 



-d/2-l 



X 



The second term in (8.1) shall subtract exactly the pole term of (8.4) divided by d, not 
the pole of (8.4), which contains a factor d = dc + 0{e). To verify this, remark that 



L 



2-D 



{bb-^ - 



cc 



-D 



=fixed, b,c 2 

which is proven by partial integration. This yields: 



,2-D 



.5) 



2-D 



b,c,g 



2^s-dc/2-l 



b,c,g 



b,c,g 



-dc/2-1 



(8.6) 



where in the last step we used the invariance by conformal mapping. This is equivalent 
to 

^' - - ' (8.7) 



which proves the desired result. We therefore can write: 




X [ln(a2'' + g^") -{2-D) ln(max(a, b, c, g))] + 0{e) (8.8) 

The method to numerically integrate (8.8) is the same as in section 7.7. We give the 
results in figure 8.1. 

In the limit D — > 1 (8.8) reduces to 



+ 9)' 



+ 2 dg 



ln(2) - 1 = -0.306852819 (8.9) 
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Figure 8.1: Numerical results for equation (8.8). The statistical error is 10 ^. 

9 Renormalization of the Wavefunction, Main Con- 
tribution 

9.1 Derivation of an analytic expression 

In this section, the calculation of ^i, defined in equation (4.31), will be discussed. We 
write: 



7(L) = = 




(9.1) 



Let us discuss the MOPE-coefficients involved. Explicit expressions are given later. The 
first is 

(9.2) has two different types of sub divergences, which are subtracted by the two counter- 
terms in (9.1). Let us symbolically write down the factorization of the MOPE-coefficients. 
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Figure 9.1: The distances in (9.1) 



If one of the dipoles is contracted first, the MOPE-coefficient factorizes as: 




^ +... (9.3) 



The first term is a relevant counterterm, which we did not mention exphcitly in (9.1). 
The second subtracts the marginal subdivergence. Note that we have again to take care 
of the tensorial structure of the factorization. 

If the two dipoles are contracted to a single dipole first, the factorization is: 



^ + 



(9.4) 



Note that any of these contractions is obtained with a combinatorial factor 2. 

Wc now give a list of the MOPE-coefficients together with the ©-functions, which 
restrict each counterterm to the sector in which the divergence appears and which by 
integration delivers the factor 1/2 from the nested integration, cf. section 7.1. First of all 



GO 



(9.5) 



with 



£2u 2 I 2i/ ^2 I / 2 I L 

/ a + a / +-(e +6 



2 c2 



<f){c 



2v _ g2i/ 



X 



X 



2y p2v 



J2u 



{c" — e 



2v 



-d/2-1 



(9.6) 



The notation for the distances follows figure 9.1. The two relevant counterterms which 
appear, when a or / is contracted first and which can symbolically be written as 



B 



(9.7) 
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are 



B 



The marginal counterterm 



(9.8) 





= Ci, C2, Ci or ^2 



(9.9) 



appears when either / is contracted first 



SD 

d + 2 



64:D 



2-D 
8D 
d + 2 



{b^ + - - d^) (if + b'- d^)h-'' - if + c'- e2)c-^) a-^('^+2) / 
(2 - Df {{f + b^- d^)b-'' - if + - e2)c-^)' a^-'^(''+2)/-^('^+^) 



u(d+2) f-v{d+2) 



X 0(/ < max:(a, 6, c)) 



X 

(9.10) 

(6^ + - - d^) ((f + - c2)e-^ - {f + d^ - 6^)^-^) a-(<^+2)_^-.(dH-2) 
64L» ^2 ~ (^-^^ + ~ ~ ^-^^ + ~ &')ci-^)' a^-(rf+2) 



X ©(/ < max(a, d, e)) 
or when a is contracted first: 



X 

(9.11) 



(6^ + - - d^) ((a^ + 6^ - c2)6-^ - (a^ + d' - e^)^-^) a-(d+2)y-.(d+2) 
64r> ~ (^""^ + ~ ~ ^""^ + ~ e2)d-^)' /i?-(<i+2)^-.(d+2) 



8D 
d + 2 



X 0(a < max(6, d, /)) 



X 

(9.12) 



2-L» 
8D 

d + 2 
"640 



^a^ + c'-b')c-^) a-'^('^+2)j-(^+2) 



(6^ + - - d^) ((^2 + e2 - d2)e-^ 

(2 - [{a' + e^- d^)e-'' - (a^ + - 52)^-0^2 ^D-.(<i+2)^-.(<i+2) 
X 6(a < max(c, e, /)) 



X 

(9.13) 



Two equivalent versions are given as one e.g. may put the counterterm for / — > on either 
endpoint of the distance /. 

The last class of countcrterms, equation (9.4), appears when either (c, d) ^ or 
(6, e) — > 0. For the first contraction, the two equivalent versions are: 



1 



Di = a 

2D 



D-ud 



c^' + d^") ''^'e(d< a)0(c< a) 



(9.14) 
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For the second contraction they are: 

D2 = -T^a''-'' (e^^ + 62-)"'^' e(6 < a)e(e < a) 

D2 = -^f"'"' i^'" + b"')''^' e(6 < /)e(e < /) (9.15) 

We now define the integrand which has to be taken in (9.1) 

F{a, b, c,d,e,f):=A-B-^ (d + Ci + C2 + C2 + Di + Di + D2 + D2) (9.16) 

We used the symmetric version of the counterterms, but could have also taken ACi + 
(1 — A)^! instead of |(Ci + Ci). We will use this freedom later in order to simplify the 
calculations. As usual this expression has to be integrated over all distances restricted to 
be smaller than L. Applying Ld/dL and mapping onto c = L = 1 results in 

J{L) = L—I{L) = [ F{a, b, c, d, e, f) max(a, b, c, d, e, f)'^' (9.17) 

uLi J c=l; a,b,d,e,f 

The integral J{L) is convergent as can be seen from a somehow tedious generalized Taylor 
expansion for the domains of possible divergences (a — ^ 0; / — 0; (c, d) 0, or as c is 
fixed for all the other distances to 00; (6, e) —>■ 0). So the term max(a, b, c, d, e, f)'"^^ can 
be dropped and the limit £ — > 0, d — > can be performed. 

We furthermore have checked that the subtracted terms Ci, Ci, C2, C2, -Di and D2 
are up to subdominant contributions in e equivalent to the terms given in equation (9.1), 
i.e. that the changes in the domain of integration do not matter (the so-called "two-loop 
miracle"). This seems to be familiar by now, but nevertheless has to be checked. 



9.2 Analytic continuation of the measure 

Now F{a, b, c, d, e, /) will be integrated numerically. The main problem is again, that the 
measure already used in section 7 and defined in (7.19) has to be improved by partial 
integration. As in subsection 7.2 we have to calculate the partial derivative applied to 
A, B, Ci, Ci, C2, 6*2, Di and D2, i.e. R, given by equation (7.22) applied to these 
terms. Hereby R will also act on the ©-functions yielding terms proportional to the 5- 
distribution. We will first discuss some ingenious methods to analytically calculate these 
"non-diagonal" terms, which work well for Di, D2 and Ci, but fail for C2 or equivalently 
€2- We then discuss another method to analytically continue the integral, which cannot 
be used for the numerics neither. 

The problem is finally solved by brute force: As in section 7 we change the prescription 
for the sectors so that the corresponding G-functions commute with R. The corrections 
are calculated numerically. 

To analyze the problem, we first write down R applied to each term. We study es- 
pecially the contributions proportional to the ^-distributions, which cannot be calculated 
numerically. In order to simplify these formulas, we use the identity vdc — 2D. 



YIA = — \- 
2D 



(2 - D)a^r'' + 2a'^ + (e^ + 6^ - - d^) (d"^ - e"^ 



X 
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-dc/2-1 



D + 2 
2D 



p-a^ + a^-p + l(e2 + 62 _ ^2 _ rf2) (^^2. _ ^2. _ ^2. ^ ^2.^ 
a2-/-^ - 1 (c2- - e^- - b^- + (f") (d-^ - e-^) 

^2-/2. _ 1 |^^2. _ g2. _ ^2. ^ ^2.y 



X 



-dc/2-2 



(9.18) 
(9.19) 



^Sn^^^ + - - (^(/2 + ^2 _ ^2)^-i. _ (^2 ^ ^2 _ g2^^-i.^ ^ 



8D 

^2 



((f + h'- d')b-^ - if + c'- e')c-^y a-'f- 



D-4 



32D 

X ©(/ < max(a, 6, c)) 



X 



2-D 

8D 



{b' + 6^-0^- d?) (if + b^- d^)b-'' - if + - e2)c-^) a-^-V-""- 



^^n^^^ ((/' + - - if + c'- e^)c-^)\-V-^-^ 



32L) 



X - [5{f - a)e(/ > b)e{f > c)+ 



X 



5(/ - 6)e(/ > a)e(/ > c) + s{f - c)e(/ > a)e(/ > &)] (9.20) 

X (l>(/2 + - c2)e-^-2 - D{f + _ 62)^-^-2 _ 46"^ + 4^"^) 

(2-L>)(2 + D) .2 + e2 _ ^2 _ ^^2 ^ g2 _ ^2)g-D _ ^^2 ^ ^2 _ ^2)^-D\ ^ 

oD v / V / 

xa-^-V-""-' 
(2-D)(2 + D) 



16D 



((/^ + - c2)e-^ - [f + d^- b^)d-'') x 
X (D(f + e^- c2)e-^-2 - D(f + d^ - b^)d-''-^ - 46"^ + 4(/-^) 

(2 ~ -D)(2 + D)2 / 2 2 2\ -D r f2 , i2 u2\j-D\'^ -2 p 

— [if +e-c)e -{f +d -b)d ) a f 



-D-4: 



32D 

X ©(/ < max(a, c?, e)) 

^ + e'-c'- d') {if + e'- c')e-^ - (f + d' - b')d-^) a-^'V" 

(2 -D)i2 + D) ^^^2 ^ g2 _ ^2y-D _ ^f2 ^ ^2 _ ^2)^-i5^2 ^-2^-i5-2 



D-2 



32D 

xj5ia-f)e{f>d)e{f>e) 



X 



(9.21) 
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2-D 



{2-D){2 + D) 





8D 






(2 


-D){2 + D) 




16 




x{a'' + d^ 


,(2 


-D){2 + D) 





{b' + - - d^) ((a^ + 6^ - 0^)6-^ - (a^ + - e^)^-^) x 



X 



X 



xO(a < max(6, ci, /)) 



+ 



2-D 
8D 



(6^ + - - d^) ((a^ + 6^ - c^)^^ - (a^ + - e^)^-^) a-^-^-^" 



X 



X 



-S(a - d)e(a > b)Q(a > /) + ]:6(a - /)e(a > b)e(a > d) 
d J 



(9.22) 



RC2 



2 -D 



V + - - («2 ^ g2 _ ^2^^-D-2^-D-2f-D 

{2-D){2 + D) 





8D 








-L>)(2 + L') 




16 




x(a^ + e^ 




-D)(2 + D) 





(6^ + - - d"") ((a^ + - ^2)6"^ - (a^ + - 6^)0-^) 

-D-4 

,2 + e2 _ _ ^^2 ^ ^2 _ ^2)^-D^ a-^'V"' X 



X 



x9(a < max(c, e, /)) 



+ 



2-D 
8D 



X 



(6^ + - - ((^2 + _ rf2)g-D _ (^2 ^ ^2 _ ^2)^-D^ ^_z^_2^-D-2 



32D 
'1 



X 



X 



1 



J{a - e)Q{a > c)e(a > /) + j6{a - f)e{a > c)e(a > e) 



RL>i = a-^(c''^ + d2'^)-'^=/'-^d-°e(d < a)e(c < a) 

+ + d^'^y^'/^ldid - a)0(c < a) 



2D 



d 



RDi = Q/-^-'(c2'^ + (i2-)-dc/2-i ^ /-^(c2- + d^'^y^'^'^-^d-''^ e(c < /)e(d < /) 



(9.23) 



(9.24) 



(9.25) 
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RD2 = a-^(e''^ + 62-)-<ic/2-ig-i3e(6 < a)Q{e < a) 

+;^«"''(e'" + b^n-'^^/^-Sie - a)Q(b < a) 
2D e 

-2^/-''(e'^ + 6'^)-'^^/'y<5(e - a)e{b < a) 



(9.26) 



First we show that in this parametrization the terms proportional to the ^-distribution in 
(9.24) and (9.26) can be calculated analytically. We demonstrate that for (9.24). 

Wc use the standard measure of equation (7.19). The vectors over which the integra- 
tion will be performed are a and d. The integration over ^3 has been improved by partial 
integration. Let us apply that to the last term in (9.24), i.e. to 

^a-^(c2^ + d^^y^/^^did - a)e(c < a) . (9.27) 

Zij-^ (Jj 

The integration over a is standard and can easily be performed, resulting in: 

— * 

The former integral over d — {x, y, z,0, . . .) was 

Sd-2 



s 



D 



/oo /• 00 /• 00 

dx dy dzz^-^F{a,b,c,d,e,f) 
-00 J —00 Jo 



and is changed to 



1 Sn-2 



2~D S 



D 



/oo /"OO rOO 

dx dy dzz^-^RF{a,b,c,d,e,f) 
-00 J —00 Jo 



(9.29) 



(9.30) 



Reversing the derivation of the measure we recognize the invariant measure in D -|- 2 
dimensions: 

1 Sn-2Sn^2 f ^_]_ r ^3 3^^ 



2-D 



02 



This can be summarized in the following formula, valid as long as the integral over d 
factorizes from the integral over a, b and c: 



/ F(a,b,c,d,e, f) — —— f RF(a,b,c,d,e, f) 



The final result is given by the integral 



2D2 



dd 
T 



-dc/2 



4L>2(2 - D 



-B 



D 



D 



2-D' 2-D 



(9.32) 



(9.33) 



where B is the standard Beta-function. It is worth to mention that this is equivalent to 



1 1 
D 2 



(9.34) 
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The factor 1/D is due to the measure, the factor 1/2 due to the fact that the R-operation 
is acting on half of the bounded distances only. The concerned reader will be able to 
reconstruct the factor 1/D from the identity 

The analogous term in (9.26) gives the same result. It nevertheless is somehow harder to 
calculate. As a first step a mapping has to be performed, which exchanges h with c and d 
with e. This mapping is possible, although the measure is not invariant, as this mapping 
does not affect the height of the tetrahedron. 

Equivalently it is possible to calculate the diagonal term of RCi, equation (9. 20). We 
find: 




(9.36) 



This time no factor 1/2 appears, since R was acting on all the distances. 

The non-diagonal part of C2 or C2 cannot be calculated by these methods. (The 
concerned reader is encouraged to try this himself.) Nor is it possible to arrange the 
counterterms in such a way that the non-diagonal terms cancel. (Warning: The measure 
is not invariant by the conformal mapping, as the partial integration by /s has been 
performed.) 

Regarding the problems with the R-operation one may ask, whether it is not better to 
use an alternative prescription for the analytic continuation. One might think of replacing 
the crucial integral 

ldfsf^-'F{a,b,c,d,e,f) (9.37) 

by 

/ d/3/3^-3F(a, b, c, d, e, f) - ffH^, b, c, dp, e^, fp) (9.38) 

where dp, Cp and fp are the projections of d, e and / onto the plane spanned by fi and 
/2 in the standard paramctrization of the measure in equation (7.19). Here the problem 
arizes that even if the length of / = (/i, /2, /s) is large, the length of fp = (/i, /2, 0) may 
tend to 0. This is a new divergence, which is integrable, but not handable by the standard 
measure in the numerical integration procedure. Thus this method cannot be used. 

9.3 Numerical integration 

We finally solve the problem by modifying as in section 7 the ©-functions and by calcu- 
lating the corrections later. We integrate RF(a, b, c, d, e, /), where 

F = A- B -C1-C2- Di- D2 (9.39) 

but with modified prescriptions for the sectors: f < d for Ci, a < c for C2, c < a for Di 
and 6 < a for L>2. We therefore have to integrate numerically: 
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Figure 9.2: Numerical results for (9.42). The first error is the statistical error, the second the correction 
for the systematic error as discussed in the text. The latter is only given, when it is not negligible. 
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32D 

X e(/ < 0?) 
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2-D 



b'^ + e^-c^- d^) (a^ + - d^jg— ^a-^'V" 



-Z)-2„-D-2 A'-D-2 



(2-D)(2 + D) 



+— 


8D ' 




X a-^-V-"" 
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16 ' 
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-L>)(2 + L') 
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2 2 

e — c 



d^) ((a^ + - d2)e-^ - (a^ + - 6^)0-^) 



X 



((^2 + g2 _ ^2^g-i? _ ^^2 ^ ^2 _ ^2^^-D^ a-^-2/-2 ^ 



j2\-D-2 



16D ' 

x0(a < c) 



X 



(9.40) 



The variable transformations are the same as discussed in section 7.3 and 7.4. The 
numerical calculations arc very difficult. Due to the complexity of the integrand, the 
numerical errors induced by the limited precision of the workstation became important 
too and limited the reduction of the statistical error. The total CPU-time was about 1000 
hours on a workstation. We obtained the numerical results summarized in figure 9.2. 

The systematic error was corrected using table 7.5. The error-bars represent the 
statistical error of the AMC-intcgration, which could not be reduced due to the lack of 
performance of the work-station, both in speed and in precision. The numerical results 
are in agreement with the analytically calculated value — | for D — > 1, discussed later. 

9.4 The correction for the first marginal counterterm 

We arranged the countcrtcrms Ci and C2 so that they both have the same correction. 
This correction is up to terms of 0{s) (for the notation cf. figure 9.4): 



c=l; a,b,d,e,f 
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2D 



if {bfb''' - cfc 



^T„-D\ ^-!/(d+2) j-u{d+2) 
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Figure 9.3: Numerical results for (9.42). The relative statistical error is 10 



+ ^-(2 - Df Cbfh-'' c/c-^)' a^-^^'^^^) 



X 



2D2 



x[e(/<c)-e(/<max(a,6,c))] 



X - [c^ — max(a, b, cY] 

£ 



X 



(9.41) 




Figure 9.4: The orientation of the vectors a, 6, c and /. 
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Expanding - [c^ — max(a, b, c)^] = ln(c) — ln(max(a, b, c)) + 0{e) and using the symmetry 
of (9.41) yields: 



Ux.V^^ ^ -cc ) a 

X [ln(c) + ln(6) - 2 ln(max(a, 6, c))] + C(e) 



X 



D-2 



fif?/^^ + c'-''' + (a' - - c2)r^c-^) 

1 n ^ ' 



a ^ X 



X [ln(c) + ln(6) - 2 ln(max(a, 6, c))] + 0(e) (9.42) 

This integral is calculated numerically using the measure given by (7.52) to (7.57) in 
section 7.7. We find the results in figure 9.3. 
For D ^ 1 an analytical calculation gives: 

i/^'d61n(6)-^^'d6 1n(6) = l (9.43) 
This is in agreement with the numerical data in figure 9.3. 

9.5 The correction for the second marginal counterterm 

Up to order e for D\ and D2 each, the following correction has to be added: 

\- [ a^-^'^ (c^" + d^")'"^' e(c < a){e{d < a) - 1) 

2D Jc=l; a,b,d,e,f ^ ' 

2D Jc=i;d>c^ ^ ' 

For — > 1 an analytical calculation gives: 

- C -^^^^ = -1 In 2 = -0.3465735903 (9.45) 

2io (l + a;)2 2 ^ ' 

The numerical data are given in figure 9.5. Numerical and analytical data fit nicely 
together. 

9.6 Analytical calculation for D 1 

We calculate now J(L = 1), equation (9.17), for D=l and 6 = 0. We give the integrals and 
the numerical results for the different regions. The combinatorical factor 4 for every pair of 
integrals compensates against the factor 1/4 from the measure. The symmetrized version 
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Figure 9.5: Numerical results for (9.44). The relative statistical error is 10" 



of the counterterms has to be used. Otherwise the 4 integrals in one equivalence class may 
not coincide. So we have to integrate ^4 — S — |(Ci + Ci + C2 + 6*2 + -Di + -Di + D2 + -D2) . 



'—{fa' + af - 2af){af - 1)-^ + lia'T' + a'T') 

+^a-2/-ie(a - /) + ^a-\a + f- 2)-'e(a - / + 1) 

+ ^a-i/-20(/ - a) + i/-i(a + / - 2)-2e(/ - a + 1) 



2 ^ ' 4: 

= / da d/0 
Jo Jo 

= 

= f da df [-lar\a - f)'' + \ {a'' + r'a'') 



(9.46) 



(9.47) 
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B = 





C =^ = 




Figure 9.6: The 3 topological inequivalent classes A, B and C and the corresponding Feynman-diagrams 
of scalar t^* theory. 



+^r\a - f)-'e{f - l)e{2f + 1 - a) 



(9.48) 



Together this gives: 



+ 



+ - 



-ln(2) + - = -0.4431... (9.49) 



9.7 The complete diagram 

In figure 9.7 we show the complete resuhs for J(l) , equation (9.17). We can verify 



£ = 



that the hmit D ^ 1 is correctly reproduced within the error-bars. 



10 Second Contribution for the Renormalization of 
the Wave Function 

The second contribution to the renormalization of the wave-function is, see (4.32): 

1 



= 



/ 




-^j +u{d + 2){ 
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Figure 9.7: Numerical results for J(l), equation (9.17). Only the statistical error is given. 



E3 



In order to determine J^2, we calculate: 




We recall that 



(10.1) 




D-v{d+2) 



Ja,b,c<L ^ ' 

/ (abb-'' - acc-^) a-^^''+^^ , (10.2) 

Ja,b,c<L ^ ' 

(a' - D {dcf c-^ (10.3) 



2D 2 
2-D 



D 
2-D 

D Ja,c<L 



-u{d+2) 



and 



u{d + 2) 



2D Ja<L 




,D-ud 



(10.4) 



Figure 10.1: The orientation of the vectors a, b and c and the position of the point x. 
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First of all one easily convinces oneself that the integral equals ior D — 1. In the 
following we give two independent derivations, which were both integrated numerically 
and which were found to coincide. This gives an additional check that the global prefactor, 
which cannot be checked analytically, is correct. 

We start with the first derivation: In order to subtract (10.4) from the first term in 
(10.2) we use the fact that 



2-D 



J a=fixed; b,c 



Therefore the first contribution is: 
v{d + 2)2-D 



2D 



Ja,b,c<L Ja<Li ^ ' 



-i/(d+2) 



(10.5) 



10.6) 



Applying \^-§i and mapping onto a — \ yields: 



u{d + 2)2-D 



Ja=l;b,c 



2D 2 

In the limit £ — > this expression simplifies to: 
{2-D){2 + D 



-'^{d+2) 1 b, c)-' - a-') (10.7) 



AD 



j (bb-^ - cc-^fa-^-^ (ln(max(a, b, c)) - ln(a)) + 0{e) (10.8) 



The second contribution is: 
2-D 



[ (abb-'' - acc-^) a-"^"-^'^ + ^— ^ / (a'-D (ac)' C') c-"" a""^^^^^ 



D 

2-D 



D 

2-D 



+ 



D 



I [(a6r^ - acc"^) - (a" - D (ac)' c'^) c 

Ja,b,c<L ^ ^ ' 

/ -/ ](a^-D{dc)'c-^)c-''a-'' 

Ja,c<L Ja,b,c<L\ ^ ' 



-D 



{d+2) 



;io.9) 



The last term becomes after application of and mapping onto a—V. 



2-D 
D 



I fa^ -L>(ac)'c-')c-°a-'^('^+')- fmax(a,c)-"-max(a,6,c)-^) (10.10) 

Ja=l; b,c^ ' e ^ ' 



In the limit £ — > it reduces to 
2-D 



/ fa2-D(ac)^c-2)c-^a-^-2fln(max(a,6,c))-ln(max(a,c))) +C(e) 

D Ja=l;b,c^ ^ ^ ' 

(10.11) 

The other term in (10.9) is: 



2-D 



D 



(abb-'^ - acc-'') -(a'-D {ac)^ c"') c 



D 



(10.12) 
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By power-counting, a pole in - can appear for a — > 0. In this limit however the terms in 
the square-brackets cancel and make the integral finite. Let us apply L-^ and map onto 
a = 1: 



2-D 



D 



J ^ ^ [[abb~^ - acc"^) - (a^ - D (ac)^ c"^) c"^] a'^'^a^ max(a, b, c 



This integral is a function of e, which has the form: 

ci£ + 0{e^) 



;io.i3) 



(10.14) 



Especially it vanishes for £ = 0. The integral thus does not change if we subtract its 
value for £ = 0. Equation (10.12) then becomes - remember that a factor ^ has to be 
reintroduced: 



2-D 
D 



abb-^ - acc-^) -(a^-D {acf c'^) c 



Ja=\; b,c '- ^ 

X - (a^ max(a, b, c)~^ — l) 



D 



X 



(10.15) 



In the limit £ ^ this simplifies to: 

^ '^abb^^ — acc^^ 



X (ln(max(a, 6, c)) - ln(a)) -\- 0{e) 
The final result is the sum of (10.8), (10.11) and (10.16): 



X 



;io.i6) 




^) +v{d + 2){ 



{2-D){2 + D) 
AD 



Ja=l; b,c ^ ' 

X ^ln(max(a, 6, c)) — ln(a)^ 



X 



+ 



2-D 



/ 



2D Ja=l;b,c 
,2 



D/A (a^ + - b^y c-') c-^(ln(a) - ln(max(a, c))) 

+ (a^ -D/A (a^ + b^- c^f b-^^ 6-^(ln(a) - ln(max(a, 6))) 

+ ((a' + b^- ^)b-'^ + (a' + - 6')c-^) (ln(max(a, 6, c)) - ln(a)) x 

xa-^-2 + 0(£) (10.17) 

Another way to treat the second term in (10.2) is to perform a partial integration: 
^ / / e(L - b)Q{L - c) oV, (6^-^ - c^-^) a-'^('^+2) 
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Figure 10.2: Numerical results for (10.17). The relative statistical error is 10 

= / / czV. [e(L - b)e{L - c)] (62-^ - c2-^) a-'^('^+2) 
= -i / / « (^(^ - b)bb-^e{L - c) + e(L - b)S{L - cjcc 

D Ja<L Jx ^ 

D Ja,c<b=L ^ ' 



1 

D Ja,b<c=L 

■-/ 

D Ja,b<c=L 



—2 

acc 



(62-^ - 



(10.18) 
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The pole in (10.3) has to be subtracted. It can be written as 

2-D 



D Ja<L:c=L 



{ac) a 



2 -v{d+2)-2-D 



and is spht into two parts 

2-D 



D 



.Ja,b<c=L Ja<c=L<b. 



[ac) a 



2 ^-u{d+2)^-2-D 



(10.19) 



;io.2o) 



which are mapped onto the same sector as in (10.18): 
2-D 



D 



f (ac)' a-'^(''+2)c-2-^ + (a6)' a-'^(''+2)&-2-^(6/c)-^e(a < b) (10.21) 

Ja,b<c=L ^ ^ 



The complete expression, i.e. (10.18) — (10.21) is: 
1 



-- f 2acc-2(62-^-c2-^)-(2-D) {ac)'c-^-''+(db)h-^-''{b/c)-'e{a<b) 

D Ja,b<c=L [ ^ ' L ^ ' 

(10.22)' 

One checks that the integral is locally convergent and that the limit e — can be 
taken. Of course this derivation is not systematic but the final result is easier to integrate 
numerically and will therefore be used in the following. It is: 



_1 / UncC-^ (62-D _ 

D Ja,b<c=L \ ^ ^ 



-(2 - D) 



(ac)' c-'-^ + (afe)' r2-^e(a < b) 



-D-2 



-I I (a' + c' - 6')c-' f 6'-^ - c^-^) + ^ 

D Ja,b<c=L \ ^ ^4 



X 



X 



(a^ + C' _ i,2)' c-^-D + (^2 ^ ^2 _ ^2^^-2-0 ^ 



a-^-^ (10.23) 



The numerical integration is performed using the measure given in equations (7.52) to 
(7.57) in section 7.7. We find the results in figure 10.2. We verify that for D — >■ 1 and for 
D ^ 2 the analytically predicted value is correctly reproduced. 



11 Contribution to the Wavefunction Renormaliza- 
tion from the 1-loop CoupHng Constant Renor- 
malization 

We will calculate the last diagram, equation (4.33), 

1 



'11. r 
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Figure 11.1: Numerical results for equation (11.5). The relative statistical error is 10 



We expand up to first order in e the following term: 

Js<LJt<L ^ ' 




-d/2 



» • 



(11.2) 



We apply L-^ and map onto s — L — 1: 



L=l 



= /~ (1 + t'^y" + r (1 + 1^')-"'' (t- - 1) 



1 r(^)r(£3) 

2-D 



2D-£ 
2-D 



+ 



(2 - DY 



Jo t 



65 



This expression can still be expanded as: 



1 r(^y 

2--Dr(S) 

t 

The final result is 



1 + £ ^ 



2D 
2-D 



- ^ 



D 



Jo t 



(2 - Dy 



2^/(2-^) ln(t) 



1 r 
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2-D 



.0 2-D 



r(^) 



2D 
2-D 



D 
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;ii.4) 



+ 



1 f' dt^D/i2-D)^^^^y2D/i2-D)^^^^^ (^^ g) 

- D)"' Jo t 



{2-D) 



The numerical results are given in figure 11.1. For D — 1, equation (11.5) is analytically 
calculated to be 

1 -ln(2) = 0.30685281944. .. (11.6) 



12 Extrapolations and Calculation of Critical Expo- 
nents 

12.1 The renormalizat ion-group functions 

As we already discussed in section 3, the renormalizat ion group /9-function and the anoma- 
lous scaling dimension v of r are in the MS-scheme that we use obtained from the variation 
of the coupling constant and the field with respect to the renormalization scale /i, keeping 
the bare couplings fixed. Written in terms of Z and they are: 



-eb 



u{b) 



bo 

2-D 



l + 6|lnZ, + f6|lnZ 



1 d_ 
2^dil 



bo 



lnZ = ^— ^- -/?(6)4lnZ 
2 2^^ 'db 



We recall the form of Z and Zb from equations (4.7), (4.26) and (4.27): 



Z = 1 + -6+ + 

£ \ e e 



= 1 



ai 



b 



'ci C2{ey 



b^ + 0{b^) 



Using equations (4.25), (5.20), (4.12) and (4.13) we obtain 



(12.1) 
(12.2) 



(12.3) 



;i2.4) 
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Figure 12.1: The functions J^i(-D) (dashed Hne) and f2(-D). The latter is given with the corresponding 
error-bars of the statistical error. 



and 



2 — D 2 — D 



(12.5) 



2 ■ AD 

For e > 0, the /^-function has a non-trivial IR-attractive fixed point for = 0. Up to 
second order in e, b* is: 



b-{e) 



1 + (2 - - - ^>f' - + 0(e^) (12.6) 



1 + ar ■ 2D(2-L')(l + ai)3 

Plugging in (12.6) in (12.5) yields: 



2 - D , 2-D ^ , 8D'a^f\ - AD(2 - D)c, + (2 - + «i) ,2 , ^(,3) (^2.7) 



+ 



2 4D(l + ai) 8D2(i + a^)3 

ci and /i were both calculated numerically as a function of D in the interval 1 < D < 2. 



12.2 The exponent u 

We can rewrite (12.7) as an s-expansion for the critical exponent u in terms of e and D. 

V* = v{D, d) = uo{D) + ui{D)s{D, d) + U2{D)s\D, d) + . . . (12.8) 

(12.9) 



2-D 

e{D,d)^2D —d 



The result for the coefficients is given in figure 12.1. In order to proceed, we shall use a 
polynomial interpolation for fi{D) and for Ci{D)/ai{D). (A polynomial interpolation for 
Ci{D) is bad as this term vanishes exponentially with l/{2 — D).) 
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Do Do 

Figure 12.2: Extrapolation of v in D and e to {D,d) = (1,3) and (2,3). The dashed line is the 1-loop 
result, the full line the 2-loop result. The exponent v is plotted as function of the expansion-point Dq. 

We now have to use these 2-loop results to calculate the critical exponent u for self- 
avoiding two-dimensional membranes (D = 2) as a function of the dimension d of space 
and to compare the results with the previous 1-loop results for u. 

As already stressed in previous works, the ^-expansion given by (12.8) cannot be used 
directly for membranes, as it can be done for polymers. There one fixes D = 1 and uses 
sophisticated resummation methods to evaluate u for d = 3 {e = 1/2) and even for d = 2 
{e = 1). Indeed, directly setting D = 2 in the ^-expansion for u gives a trivial, but absurd, 
result, since all the terms ^'^(2) of the e-expansion vanish! Moreover when = 2,6 = 4 
irrespective of the value of d. This simply means that the point D = 2, e = (which 
corresponds to tic = oo) is a singular point and that it is not possible to perform a direct s- 
expansion around it. Instead, one may perform a similar expansion, starting from another 
point {Dq, do = 2-^^o ) '^^ critical curve {e = 0). This approach has already been used 
in [30] to extrapolate the 1-loop results. It introduces two different kinds of arbitrariness 
in the extrapolations. Firstly one is a priori free to chose any starting point on the critical 




Do Do 

Figure 12.3: Extrapolation in D and d to {D,d) = (1,3) and (2,3). 
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j^:_(r>,d) = (i,i) 





(D,d) = (l,3) 



Figure 12.4: Extrapolation for polymers in D and Dc{d) to {D,d) = (1,1), (1,2), (1,3) 



curve e = (or at least any point in some subset of this curve). Secondly there is an 
arbitrariness in the "path of extrapolation" which goes from the critical curve to the 
physical point {D = 2,d). Of course if one knows an exact resummation procedure of the 
full series in e, one should obtain the same result for any starting point and for any path 
of extrapolation. However, if one has a finite number of terms of the e-expansion, even 
for an adequate resummation method, the final result will depend on this arbitrariness. 

In [30] a minimal sensitivity method was introduced by Hwa to analyze the 1-loop 
results. He chose a given extrapolation method and selected the starting point {Do,do) 
on the critical curve which gives an estimate for i/, z/(Do) which is the less sensitive to the 
choice of the starting point. This method works well for polymers and gives interesting 
results for membranes. However, when one goes to 2-loops, it delivers largely varying or 
even diverging estimates for i^(i^o) and in addition his choice for the extrapolation path 
is somewhat arbitrary. 

We shall use a generahzation of the methods introduced in [30]. Note that the ex- 
pansion (12.8) is exact in D and of order 2 in £, thus it can be expanded up to order 2 
both in D — Dq and e. Now we can change our extrapolation path through any invertible 
transformation {x, y} = {x{D,s),y{D,e)}. One can express D and e as function of x and 
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y and re-expand u up to order 2 in x and y around the point {xq, yo) on the critical curve. 

p{D, e) = u{x, y) = z>o,o(a;o, yo) + Axui^ixo, yo) + A|/z>o,i(xo, yo) 

+{Ax)'^i>2,o{xo,yo) + 2AxAyi>i,iixo,yo) + (Ay)^i>o,2(a;o, Z/o) H 

Ax = x-Xo, Ay = y-yo (12.10) 
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Figure 12.6: Extrapolation in D and l/{d + 2) to = (1,1), (1,2), (1,3), (2,2), (2,3) and (2,20) 



The goal is to find an optimal choice of variables {x,y}. Our guidelines for such a 
choice are the following: (i) the estimate for u should depend "the least" on the choice 
of the expansion point on the critical curve, (ii) it should reproduce well the known 
result for polymers {D = 1), {Hi) for membranes {D = 2) the d —> oo limit should not 
be singular and we should get results close to those obtained by a Gaussian variational 
approximation. This last point is not arbitrary and will be justified below. It turns out 
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j^:_(r>,d) = (i,i) 




(D,d) = (l,3) 



Figure 12.7: Extrapolation for polymers in Dc{d) and e to {D,d) = (1,1), (1,2) and (1,3) 



to be quite stringent. 

Finally, we must choose some resummation procedure to extrapolate u from the know- 
ledge of the series (12.10) up to order 2. Since we have only a few terms and since we do 
not have insight in the large order behavior of these series or in the analytical structure 
of the resummed series, we cannot use sophisticated resummation methods (for instance 
those based on Borel transforms). Therefore we shall always use the truncated series at 
order 2 and boldly sum its terms. 

Let us now discuss possible extrapolation variables. The simplest choice is to take D 
and e. This works well for polymers {D — 1), since both at 1- and 2-loop order we get 
results which are quite stable with respect to Dq. However this gives very poor results 
for membranes {D = 2), since both at 1- and 2-loop order the results depend very much 
on Dq. This could be expected, since in this case e = 4 independent of d. See figure 12.2. 
In the case of the membrane, no prediction is possible. 

Another possibility is to expand in D and d. This expansion gives reasonable results 
for polymers, but poor results for membranes. This can be seen by looking on the results 
of figure 12.3 and is not surprising, since if we apply this extrapolation method to the 
Gaussian variational estimate z/^ar = 2D /d it also gives poor results, although the result 
i^var is expected to be close to the exact u at large d (as argued in subsection 12.3). 
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Figure 12.8: Extrapolation in Da{d) and £ to {D,d) = (2,2), (2,3), (2,4), (2,8), (2,15) and (2,20). 



A more interesting choice is to use D and Dc{d) = 2d/ (4 + d) or equivalently D 
and + 4). This expansion has the advantage to represent the critical curve Dq, 
dc{Do) as a straight hne. For polymers in 3 dimensions this method delivers a remarkable 
broad plateau, i.e. the extrapolated value of u is relatively independent of the expansion- 
point Do- For polymers in 2 and 1 dimensions we still have a plateau when Dq 2, 
which delivers 2-loop extrapolations close to the exact results {u = 0.75 and = 1 
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Figure 12.9: Results of the numerical extrapolations for v. If not stated otherwise the error is ±1 in the 
last digit, (a;, y) indicates the expansion-parameters x and y as discussed in the text, min, max and mean 
are the minimum and the maximum of the plateau and their mean- value respectively. 

respectively) (see figure 12.4). For membranes, this method also delivers interesting results 
(sec figure 12.5). For large d, wc find a stable plateau at 2- loops when Dq — > 2, which 
gives for v a result very close to the variational estimate z/yar = 4/(i, while at smaller Dq, 
the plateau stops and the 2-loop estimate for v increases sharply as i^o — 1- For small 
d, the Dq ^ 2 plateau becomes an oscillatory region, still followed by a sharp increase 
for Dq ^ 1. In this case, to calculate u, we use the minimum and the maximum of 
v{Dq). Whereas the first is expected to be an underestimation, the second is expected 
to be an overestimation. We also give their mean-value, cf. figures 12.9 and 12.5. One 
might think of developing in D and l/(d + c), with c ^ A. c = is suggested by the 
variational ansatz. In fact we prefer to take c = 2, which is suggested by the prediction 
of the Flory argument (see subsectionl2.4). We obtain similar estimates as above for d 
large, but larger variations for smaller d. For this case we evaluated v by the request that 
the second order corrections should vanish. We find the results in figure 12.9 and 12.6. 
The predictions for u are good in the known cases and reasonable for membranes in 3 
dimensions. 

Another promising method is the expansion in e and Dc{d). This expansion is also 
regular for D ^ 2 and is perhaps more in the spirit of an e-expansion. Let us discuss the 
features of this expansion in more detail. See figure 12.8. For polymers in 3 dimensions we 
find the flattest plateau of all extrapolation methods. For membranes in 3 dimensions the 
prediction at 1-loop order (dashed line) is essentially independent of the expansion point. 
For membranes in large dimension, at 2-loop order the estimate starts from the 1-loop 
result at D = 2, grows until it reaches a plateau, where u ^ u^ar and then grows rapidly 
again. For smaller d, there is still a plateau and in order to extract from figure 12.8, 
one uses the maximum and the minimum of the plateau. Their mean is an estimate for 
v, their difference an estimate of the error in this expansion scheme. 

The results for i/ from the various extrapolations are summarized in figure 12.9. 
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12.3 Variational method and perturbation expansion 

In the various extrapolation schemes that we used, we have seen that the plateau structure 
for Do — > 2 becomes clearer when the space dimension d is large and that the correspond- 
ing estimates for 1/ are close to the value i^var — ^/d obtained from a variational ansatz. In 
order to understand this phenomenon and to have a better understanding of the plateau 
structure of the extrapolations, we shall discuss the status of the variational method. 

Using a Gaussian variational ansatz [39] , the exponent v for the crumpled phase of a 
self- avoiding D-dimensional tethered membrane is found to be 

2D 

i^var = -J- . (12.11) 

It was noticed in [31] that in the case of membranes with long range interactions the 
variational estimate for u is exact and can be reproduced easily by simply assuming 
that there is no coupling constant renormalization, i.e. that one can take — 1 in the 
renormalized Hamiltonian (3.13). This last assumption can be proven for membranes 
with long range forces [31, 32]. Let us rewrite the full dimension i/(6) in a way which 
makes this point clear. Starting from the definitions of the RG-functions f3{b) (12.1) and 
iy{b) (12.2), we reexpress ^InZ in terms of /3(6) and ^InZf,, 

Inserting this into (12.2) gives 

At the IR fixed point b = b*, the second term of the r.h.s. of (12.13) vanishes. The last 
one does not vanish in general, but vanishes if = 1. In this case, we get the variational 
result (12.11). 

Using this observation, one can understand why the 1-loop and 2-loop extrapolations 
for u coincide with z/^ar for large d. Taking do 00 on the critical line amounts to take 
Dq ^ 2. The limit d 00 corresponds thus to the limit D ^ 2 for the counterterms. 
The 1-loop wave-function counterterm is given by the residue 




-^j ^1 as D^2 (12.14) 

while the 1-loop coupling constant counterterm is given by the residue 

^2-2^/(2-D) as ^^2 (12.15) 



which is exponentially smaller than (12.14) when D — > 2. A similar exponential factor ap- 
pears for the 2-loop coupling constant counterterm compared to the 2-loop wave-function 
counterterm (this can be checked from the analytical expressions and the numerical re- 
sults) . 
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Figure 12.10: Extrapolation of i/yar = 2£)/d for membranes in {£)c(d),e} (left), and {£'c(-D),£'} (right), 
to (£',c?) = (2,3), (2,10) and (2,20). On the first figure the labels 1-4 indicate the order of the e- 
expansion i.e. O(e^) - O(e^). On all figures, the full line is the result at order the other lines are the 
results at order £, and 



When looking at the general structure of the divergent diagrams at iV-loop order, we 
can argue that this phenomenon will persist, but we have no rigorous proof. However, 
if this exponential bound InZ;, ^ InZ as — > 2 is correct, this means that v — v^^x ~ 
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Figure 12.11: Results for v using the numerical extrapolations for vd. If not stated otherwise the error 
is ±1 in the last digit. 



exp(— cst./d) when d — > oo. 

Assuming the variational estimate (12.11) to be a good approximation for large d, 
it is interesting to test the various extrapolation methods that we have used to get the 
2- loop results. The principle is to start from the exact formula for i/yar, to write an £-like- 
expansion around the critical curve £ = 0, to truncate it at a fixed order and to resum 
the result as done previously. The results of such a resummation, using the extrapolation 
in are presented on figure 12.10. Let us note the following points: 

1. The extrapolations at 1- and 2- loop orders of i^var are indeed very similar to the 
extrapolations of v at large d. In particular one recovers the same plateau structure 
and it should be noted that already the 2-loop optimal estimate for z/yar, obtained 
by the minimal sensitivity method, gives the exact VysJ. This is a strong point for 
this method, when applied for smaller d. 

2. One should note that the ^-expansion for i/yar is convergent, but that is has a finite 
radius of convergence. As a consequence, one can show that the A'"-loop extrapola- 
tions for i/yar converge towards the exact result when — > oo only in a finite range 
of starting points Dq on the critical curve. This range is explicitly 

< Do < 2 (12.16) 



8 + (i 



3. The large variation of the 2-loop estimates for v a& becomes small reflects the 
fact that we are outside of the range of convergence of the ^-expansion. This is 
clear when one compares the 2-loop estimate with the 3- and 4-loop estimates on 
figure 12.10. We expect that this is still true for smaller d and that the optimal 
values for Dq are still in the domain of confidence of the ^-expansion. 

We now present a new e-cxpansion, well suited to large d, which is suggested by the 
expression (12.13) for vih). Indeed, when evaluating the exponent v* — v{h*^ by (12.13), 
we see that what we really expand in e is not v but vd: 

v*d ^2D+ ('/?(6)^lnZ,(6) 

^ ^ 6=6* 

= 2L> + 0(£) (12.17) 
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Figure 12.12: Extrapolation of ud in D and Dc{d) to {D, d) = (1,2), (1,3), (2,2), (2,3), (2,6), (2,20). 



Thus we may perform the ^-expansion for ud rather than for u. This new expansion has 
the advantage that it starts at order £° from the result predicted by the variational ansatz. 
The same extrapolation methods used for v can be used for vd. 

The results are given in figures 12.11 and 12.12. Wc find that for polymers {D — 1), 
when using the {D,Dc{d)) variables, the minimal sensitivity method gives poor results, 
but that the criterion to take as optimal Dq the point where the second order correction 
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vanishes gives good results. We use the same criterion for membranes {D — 2). 



12.4 Expansion around Flory's estimation 

In the last section we have seen that the e-expansion for ud is in fact an ^-expansion 
around the variational ansatz, i/d — 2D + 0{e). Another stimulating idea is to perform 
a similar expansion around the prediction made by Flory's argument. It is well known 

that the Flory result for polymers z/piory = 3/(2 + (i) is simply obtained by assuming that 
the elastic term and the contact interaction term in the Edwards Hamiltonian scale in 
the same way with the internal size of the polymer. The same scaling assumption for 
membranes leads to the prediction 

2 + D , , 

It is possible to perform an ^-expansion around i/piory by simply expanding ^{d + 2). 
Indeed, we can set 



^Fiory = y Z Zh , = yZ/Zb (12.19) 

and eliminate Zpiory from the system of equations (12.1) and (12.2). We obtain 



^*{d + 2) = D + 2+(m^^i^{^ 



(12.20) 

b=b* 



This makes clear that if the wave-function and coupling-constant renormalizations are 
the same (more precisely if Zf^/Z stays finite at the IR fixed point b*) the Flory result 
becomes exact. Moreover, the £-expansion of i/(d -|- 2) is clearly an £-expansion around 

i^Flory 

This expansion seems to be the most satisfying numerically. In particular, the method 
of minimal sensitivity and that of minimizing the second order term give generally close 
results. We therefore give in the following plots for all interesting combinations of vari- 
ables and dimensions {D, d). Let us stress that good expansion-parameters for i/ are not 
necessarily good for ^{d + 2) and vice versa. E.g. the expansion in D and d is bad for 
u but works quite well (although not optimal) for ^{d + 2). We study the expansion in 
{D, d), {d, e), {D, Dc{d)) and {Dc{d),e). The results of the extrapolations are collected in 
figure 12.13. 
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Figure 12.13: Results for i/ from the numerical extrapolations for i'{d + 2) at second order, "min", "max" 
and "mean" denote respectively the estimate at the minimum of the plateau, at the maximum of the 
plateau and their mean value, "left crossing" , "right crossing" and "mean crossing" denote respectively 
the estimate at the leftmost point where the second order correction vanishes, at the rightmost point and 
their mean value. If not stated otherwise the error is ±1 in the last digit. 
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12.5 Summary of the 2-loop extrapolations for v. 



Let us summarize. In figure 12.18 we represent the results of a 2-loop extrapolation for v 
in the case of membranes {D — 2) in d dimensions {2 < d < 20). We see that for d — > oo 




Figure 12.18: Extrapolation of the 2-loop results in d and e for membranes D — 2 in d dimensions, 
using (12.20) (squares). The sohd line is the prediction made by Flory's theory, the dashed hne by the 
variational ansatz. 

the prediction of the Gaussian variational method becomes exact, as argued above. For 
small d, the prediction made by Flory's argument is close to our results. This is a non- 
trivial statement, since the membrane case corresponds to £ = 4 and in comparison with 
polymers in d = 3, where e = 1/2, the 2-loop corrections were expected to be large. In 
fact we have found that the 2-loop corrections arc small when one expands around the 
critical curve £ = for an adequate range of D ^ 1.5 (depending slightly on d and on 
the choice of variables) and a suitable choice of extrapolation variables. In this case the 
2-loop corrections are even smaller than the 1-loop corrections and allow for more reliable 
extrapolations to £ = 4. 



12.6 Other critical exponents 



The 2-loop calculations presented in this paper allow in principle to compute other scal- 
ing exponents for self-avoiding tethered membranes. The first exponent is the so-called 
correction to scaling exponent u which governs the corrections to the large L scaling be- 
haviour. It is known that this exponent is given by the slope of the /3-function at the IR 
fixed point b* 

oj = ^J3{h) (12.21) 



d¥ 



b=b* 
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Its £-expansion is given by 



i2-DMD)+2DMD) ^ (12.22) 

(For the definition of ci{D) and fi{D) cp. (4.34) and (4.28)). 

We have tried to use the extrapolation methods that we have developed for u to get 
2- loop estimates for uj. Unfortunately, it turns out that the term of order is always 
very large compared to the term of order e and there is no domain along the critical curve 
where a reliable estimate can be extracted from the £-expansion. Thus the situation for 
the exponent cu seems to be very different from that for u. 

Another scaling exponent which can be obtained from our calculations is the so-called 
bulk contact exponent 62- For a general introduction to contact exponents for polymers 
and membranes we refer to [41]. A detailed discussion of the calculation of contact 
exponents within the Edwards model for self-avoiding membranes will be given in [32] . We 
simply recall the basic results here. The contact exponent 62 is related to the probabihty 
to find two fixed points Xi and X2 inside the membrane at a relative distance r — |f| in 
external d-dimensional space: 

P(r; X,, X2) = (Hf- (f(xi) - f{x2)))) (12.23) 

For a large membrane, P is expected to take the scaling form 

P(r; xi, X2) = R^^F{r/R^2) (12.24) 

where R12 is the mean distance between xi and X2. 

R\2 = \{{r{x,)-r{x2)f) (12.25) 
The contact exponent 62 is given by the small r behavior of the scaling function F 

when r ^ (12.26) 



M12J \Rl2. 

02 is simply related to the scaling dimension UJ12 of the 2-membrane contact operator 

512(^^1,0:2) = 5"(fi(a:i)-f2(x2)) (12.27) 

in the model of two independent self-avoiding membranes. This model is described by the 
Hamiltonian 

n=(f l{vn{x^yf+f ^^k^r2{x2)r) 

+&(/ / 5\n{x,) - vM) I I ~5\r2{x2) - r2{y2))) 

\JxieMi JyieM-i JX2&M2 ^2/2GM2 / 



+ 2t / 6''{n{x,)-f2{y2)) (12.28) 
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This model can be made UV-finite at £ = by the same renormahzation factors for 
r and h as the 1-membrane model, i.e. by replacing in (12.28) r — > Z{by^^f and b — > 
bZh{b)Z{b)'^^^ ji^ , but with an additional renormahzation for the inter-membrane coupling 
t tZt{b^t)Z{bY^'^ jjL^ . The new counter-term Zt can be calculated in terms of the same 
divergent diagrams as those which contribute to Zf,, but with different numerical factors. 
In particular, one can show that when t — b, Zt{b,t = b) = Zb{b), so that the symmetric 
2-membranes model reduces to the 1-membrane model. 

As a consequence of this formalism, one can define a new RG function, fit 



(3t{b,t) = i^^t 



(12.29) 

bare b, t 
fixed 



calculate the RG-fiow in the {b,t) plane and check that {b,t) = {b*,b*) is the IR stable 
fixed point which governs the scaling behavior of a large membrane. It turns out that 
the $2 contact exponent and the anomalous dimension 0012 of the contact operator 5i2 are 
related to the ^-derivative of the Pt function at the IR-fixed point by 



d 
ot 



(12.30) 



As a result, the e-expansion for 6*2 involves the same diagrams as uj and but in a different 
combination. The result is somehow complicated, so let us simply write the counterterm 
Zi and compare it with (already given in sections 3 and 4). 

z \ z j 

(12.31) 

We have tried to use the extrapolation methods that we have developed for v to get 
2- loop estimates for ^2- As for u, it turns out that the term of order £^ is always very 
large compared to the term of order z and that there is no domain along the critical curve 
where a reliable estimate can be extracted from the £-expansion. 



12.7 Scaling for membranes at the ^-point (tri-critical point) 

In [34] it was shown that the scaling behavior of membranes at the tri-critical ^-point 
is dominated by the modified 2-point interaction 5"{r{x) — f{y)), which is repulsive at 
short distance but attractive at larger distance. Membranes at the ^-point are thus 
quite different from polymers, where the 3-point repulsive interactions dominate. As a 
consequence, the scaling exponent at the ^-point has a different £-expansion, which 
can be computed analytically at first order and is found to be [32, 34] 
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with 

s' ^2>D-2- d^^-^ (12.33) 

This expansion is again correct up to first order in d or e' and exact in D. We can use the 
extrapolation methods developed in this paper to evaluate the exponent u from the 1-loop 
results. The extrapolation method proposed by Hwa [30] does not work in this case. Two 
variants of our method deliver reasonable results: The expansion in s' and Dc{d) and the 
expansion in D and Dc{d). The 1-loop estimate for membranes in 3 dimensions is (cf. 
figure 12.19): 

i/^ = 0.42 ± 0.08 (12.34) 

Due to our experiences with the former extrapolations, we expect that the 1-loop results 
will be an underestimation. They may be compared to the Flory-result 



- 1^ , (12.35) 



which evaluates to 0.57 in the case of membranes in 3 dimensions. 



12.8 Comparison with numerical simulations and experiments 

The numerical study of two-dimensional self-avoiding tethered membranes imbedded into 
three dimensions was started in 86 by Kantor, Kardar and Nelson [11, 12]. They found 
for the exponent u a value close to the prediction made by Flory's argument. These 
simulations as most of the following use a system composed of balls and springs (spring 
and bead model). Self-avoidance is effective between the balls. It was a surprise when it 
was found [14, 9] that the simulations of larger membranes obtain a flat phase. Abraham 
and Nelson[10] explained this result by suggesting that an effective bending rigidity is 
induced by the geometric constraints of the model (there is a maximal angle for which 
neighbouring faces can bend due to the finite size of the balls). If this indTiccd effective 
rigidity is larger than the critical rigidity, where the crumpling transition (transition to a 
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flat phase induced by the bending rigidity without self- avoidance) occurs, the membrane 
wiU always be in the flat phase. 

After these studies, several attempts were made to reduce the effective rigidity of 
the membrane in order to observe a crumpled phase, by using smaller balls [15] or by 
bond-dilution [22, 23]. None of these attempts was successful. Another possibility to 
reduce the rigidity is to impose self-avoidance not between balls but between faces of the 
lattice. Using this method Baumgartner et al. found a crumpled swollen phase. KroU 
and Gompper later declared that also in this model the membrane is flat [19]. Finally 
let us remark that a crumpled phase has been found by balancing the induced rigidity 
by long-range attractive interactions [16, 17]. We still want to mention the simulations of 
self-avoiding tethered membranes in a space of 4, 5, 6 and 8 dimensions. Grest [20] found 
that in 4 dimensions, the membrane is flat, but is crumpled in dimensions higher than 4. 
Subsequent simulations by Barsky and Plischke [21] conflrm this conclusion. 

To summarize: For 3 dimensions most numerical simulations indicate that the mem- 
brane is in the flat phase, but the results are still not fully conclusive. The simulations 
which flnd a crumpled phase give values for i/ which are close to our analytical estimates. 
The situation seems to be similar in 4 dimensions. 

A few experimental studies of tethered membranes have also been performed, but the 
situation is not clear neither. The best studied system is graphite oxide, i.e. a mono- 
layer of carbon-atoms. The flrst experiments by Hwa et al. [24] found a phase with 
fractal dimension df = 2.4. This result is contested in [25] by direct electron-microscopy 
methods, although the latter authors obtain within the error-bars the same data from 
their diffraction experiments. It is not clear if these membranes are sloppy enough or if 
their internal stiffness is sufficiently large to induce in itself the transition towards a flat 
phase. Another system studied is the spectrin network of red blood cells [26] , but the role 
of disorder (which may induce a wrinkling transition) seems here important. In summary, 
the experiments are not yet helpful to clarify the situation. 

12.9 A possible scenario for low dimensions 

The fact that most numerical simulations only flnd a flat phase (with u — 1) ior self- 
avoiding tethered membranes in 3 dimensions, while our analytical estimates leads to 
< 1, has still to be explained. As mentioned above, a possible explanation is that 
the analytical calculations apply to floppy membranes (with very small bending rigidity), 
while '"real systems" studied numerically or experimentally are rigid enough to cross the 
crumpling transition barrier and to stay in the flat phase. 

However, the fact that reducing the rigidity does not induce the crumpling transition 
suggest another, more drastic possibility. It is known that for phantom tethered mem- 
branes, the lower critical dimension di, below which the crumpled phase does not exist, 
lies between 1 and 2 (1 < < 2). It is possible that, when self-avoidance is taken into 
account, the lower critical dimension di for the crumpling transition moves upwards to 
about 3. In this case self-avoiding tethered membranes are always flat, however small the 
rigidity is. 

Let us give heuristic arguments that this is indeed the case in 3 dimensions. Let us 
start from phantom membranes (without self-avoidance). The fractal exponent is then 
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in the crumpled phase, 1 in the flat phase, and equal to 




(12.36) 



at the crumpling transition [42] . This last estimate is the result of a large d expansion. 
Its applicability to low dimensions is thus not clear a priori, but numerical simulations 
[8] show that even in 3 dimensions this approximation is reasonable = 2/3). Let us 
now ask whether self-avoidance is relevant at the crumpling transition. By naive power- 
counting we find that this is the case if 

V ^2-2-u^d>0 (12.37) 

i.e. with (12.36) for 

d<5 . (12.38) 

For "D > we expect that the fractal exponent Vc at the crumpling transition is different 
with or without self- avoidance. But let us assume that self-avoidance does not change 
drastically the exponent u^. at the crumpling transition, i.e. that 

^ u,+sA (12.39) 

where the second exponent is the exponent v at the crumpling transition in the presence of 
self- avoidance. This is certainly true for V small. As we expect that the radius of gyration 
of a self-avoiding membrane with rigidity is an increasing function of the bending rigidity 
modulus K, until we reach the crumpling transition, it is simple to deduce the inequality 
between the fractal exponent of self-avoiding membranes in the crumpled phase usa and 
at the crumpling transition Vc+sa 

l^c+SA < l^SA . (12.40) 

Assuming that (12.39) holds, using the estimate (12.36) for Uc+sa and our two-loop esti- 
mates for Us A, we find that this basic inequahty is violated if 

d<di^ 3.8 (12.41) 

Our interpretation is that di is nothing but the lower critical dimension of the crumpling 
transition with self-avoidance, and that for d < di the membrane is always in the flat, 
rigid phase. 

We want to emphasize that this line of arguments is still somehow speculative. It is 
however tempting to compare our crude estimate of di ~ 3.8 with the numerical simu- 
lations in 3, 4 and 5 dimensions which yield 4 < ci; < 5. To confirm or disapprove this 
scenario, one must take into account the effect of the bending rigidity in our renormaliza- 
tion group calculations. 



13 Conclusions 

In this article we presented the first renormalization group calculation at 2-loop order for 
self-avoiding fiexible tethered membranes. These second order corrections were found to 
be surprisingly small if one uses an adequate extrapolation scheme. 
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We were able to clarify the status of the Gaussian variational method and to show 
that it becomes exact for d — > oo. For low dimensions the 2- loop results are in good 
agreement with the prediction made by Flory's approximation, although systematically 
slightly larger. 

In order to improve these results, one should understand if the plateau phenomenon 
observed at 2-loop order persists to higher orders, and one should control the general large 
order behavior of perturbation theory for this model. Another important issue is whether 
the IR-fixed point studied here is stable towards perturbation by bending rigidity. Indeed, 
for small enough d this might destabilize the crumpled phase and explain why numerical 
simulations in d = 3 normally see a flat phase. We gave additional arguments which 
corroborate this scenario. 
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A Normalizations 

We use peculiar normalizations in order to simplify the calculations. First of all, we 
normalize the integration measure of the internal space as {Sd is the volume of the D- 
dimensional unit-sphere) 




2 




71- 




This provides 




(A.2) 



The ^-distribution is normalized according to 



Jp 



(A.3) 



with 




p 



(A.4) 



to have 




p 



(A.5) 



Using for the free Hamiltonian 




(A.6) 



delivers 




(A.7) 
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Figure B.l: The distances and the points in (B.l) 

B Example of the MOPE 

We give as an explicit example of the MOPE the derivation of (6.2). 



xf 'y\ xf 'j/2 xf ' y-i 

— III' o**^''*^^^) • • o*J"'(^2) • • J'l'^i^i) ■ ■ o-^kr{yi) . . ^-ipr{y2) . . ^-iqriys) . 
J k Jp Jq 



^iKr(xi) . . Qtpr(X2) . . ^iqryxi) . . ^-tHr^yi) . . ^-ipr(y2) . . ^-iqr^ys) . j-g -j^^ 

These exponentials shall be contracted hke 

• — — — • ^ ^O'^ ^^'^^ 

We therefore use the OPE for the points xi, X2 and X3, supposed the differences between 
these points become small: 

. gifer(a;i) . . ^pr{x2) . . ^qr{xz) . _ . Qikr{x\)+ipr{x2)+iqr{x3) . ^kpo?'^ +kqb'^" +pqc'^" g'^ 

The new variables for the distances between the points are given in figure B.l. An 
analogous expansion is valid for yi, y2 and yz- In order to retain only the most important 
contribution, we expand 

. ^ikr{xi)+ipr{x2)-]-iqr{xz) . _ . Qi{k+p+q)r{{xi+X2+X3)/3) ^-^ _j_ Q(J\/j'^^ ■ 4) 

and neglect the contributions of order C(Vr) because they arc proportional to irrelevant 
operators. After a shift in the integration-variable q, 

q — >q-k-p , 

equation (B.l) becomes: 

/ • Q^qr {{xi+X2+X3)/3) . . g-igr((yi+2/2+J/3)/3) . ^ 
Jo 

" ^kpia^-'+(Pn+k(g-k-p){b''-'+e^'')+Pig-k-p){c^''+f'') 



X Qkp[a--+d--)+k{q-k-p)(b-''+e-")+p(q-k-p)(c-"+f--) 
Jk Jp 

The integral over q yields the (5-distribution plus higher derivatives of this distribution. 
The latter are irrelevant operators and can be neglected. As they come from the expansion 
of the last exponential factor in q, we only have to retain the last factor, evaluated at 
q = 0. This gives: 



f f ^kpia'^-'+d?-' 
Jk Jp 



)-k{k+p){b'^''+e^'')-p{k+p){c'^''+P'') 

d/2 



(B.6) 
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This can still be factorized as is known from ancient Heron: 



-d/2 



(B.7) 



C Equation of Motion 



The equation of motion reflects the invariance of the functional integral under an change 
of variables for the fields. The expectation value of an observable O in the free theory is: 

JPHOe-^y (CI) 

We now perform a global rescaling of r{x): 

r{x) — > {1 + K)r{x) . (C.2) 

The expectation value of O, equation (C.l) is unchanged. The exphcit form becomes up 
to first order in k: 

_ JP M o (1 + « 10],.) (i - a? /. ^) ^ 

* /I'H(l-^4+)e-A/.+ 

[O]^ is the canonical dimension of the operator O, measured in units of r, i.e. that [r]^ — 1. 
Calculating the difference of (C.l) and (C.3) gives: 

/ r \ conn / y conn 

{Oj^l =.ia(0\ (C.4) 
For several operators we have: 

/ r \ conn / \ conn 

(0102 J =u{[0^l + [02l)(^0^02)^ (C.5) 

and in particular 

/\ conn / \ conn 

/ r \ conn / \ conn 

O^J^j^ =(2 + PUi^(0^)^ (C.7) 

These relations are equivalently valid for non-connected expectation values. (To prove 
this, remark that (-j-)o = 0.) 
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We can now apply these equations to calculate the pole-term of { '-^ ^ ' 

The equation of motion yields 



. . . -f ) = -2vd 

I 

On the l.h.s. the possible divergences proportional to 
integrals: 



(C.8) 

come from the following 

1:;^^ (C.9) 



^)) +2 

- / 

In the second term, -j- has to have some distance L from the endpoints of : •^^^j^ ' • 
The pole-terms proportional to • • appearing on the r.h.s. of (C.8) and in (C.9) 



are: 



+ C»(£° 

+ e'(£°) 



(C.IO) 



In the last equation, we used that 



+ 



(Cll) 



and that 



/// 



III 



^ -2 





(CT2) 



+ 0(£°) 
(C.13) 



V; ; (C.14) 





The second term has no contribution proportional to 
The first is evaluated using the equation of motion: 



• • ■ 



m • 



and thus can be neglected. 

(C.15) 



Together this yields (C.12). 

Finally using equations (C.8) to (C.12) gives: 



• • 



• • 



(C.16) 



We further have to show that 




^) . (C.17) 
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The equation of motion yields 



• • 



=-ud 





• • 



(C.18) 



The integral on the l.h.s. of (C.18) can be decomposed as: 




(C.19) 



The divergencies of the r.h.s. of (CIS) and of the first term in (C.19) are: 



// 





+ 



The divergence of the second term in (C.19) is extracted as 







(C.20) 
(C.21) 




The second factor on the r.h.s. of (C.22) is 



(C.22) 



(C.23) 



where again the equations of motion are used. Using (C.18) to (C.23) yields (C.17). 
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